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In  the  previous  contract  period,  experimental  as  well  as  theoretical  re¬ 
search  has  been  carried  out.  The  former  has  concentrated  on  the  development 
of  a  novel  testing  apparatus  which  promises  to  reveal  hitherto  inaccessible 
details  on  the  surface  condition  of  materials  during  slidinq.  It  employs  si¬ 
multaneous  measurements  of  the  coefficient  of  friction  and  the  electrical  con¬ 
tact  resistance,  using  a' "soft"  set-up  and  following  the  momentary  variations. 
This  research  has  not  yet  yielded  reportable  results  but  progress  is  most  en¬ 
couraging. 

Theoretical  work  has  been  devoted  to  theory  of  plastic  deformation,  and 
a  review  on  workhardening  theory  has  been  written  due  to  be  published  in  Trans 
Met.  Soc.  While  this  is  progressively  leading  to  a  deeper  understanding  of 
metal  plasticity,  and  in  this  much  is  valuable,  no  particularly  noteworthy  new 
results  have  emerged. 

Of  by  far  the  greater  interest  in  connection  with  the  present  contract 
theme,  namely  friction  and  wear  and  surface  processes  governing  these,  has 
been  the  research  embodied  in  the  two  extensive  papers  which  are  appended 
herewith.  As  is  seen  from  these,  it  has  been  possible  to  derive  rather  sim¬ 
ple  analytical,  closed  equations  for  the  evaluation  of  the  momentary  temper¬ 
ature  at  contact  spots  during  sliding.  This  theory  has  been  extended  to  the 
case  of  concurrent  friction  and  Joule  heating,  as  presented  in  ttre-thi rd ,•  •' 
-shorter  enclosed  paper,  which  has  been  submitted  to  the  WEAR  OF  MATERIALS  - 
1985  conference. 
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DE-MYSTIFYING  FLASH  TEMPERATURES  -  II  APPLICATION  TO  HEATING  THROUGH  FRICTION^ 

AT  PLASTICALLY  DEFORMED  CONTACT  SPOTS 

by 

D.  Kuhlmann-Wilsdorf 
Department  of  Materials  Science 
University  of  Virginia,  Charlottesville,  VA  22901 

Abstract 

In  Part  I  expressions  for  the  flash  temperatures  at  elliptical  contact 
spots,  AT^,  were  derived.  They  were  uniformly  found  to  be  proportional  to 
q,  the  rate  of  heat  evolution  per  unit  area  of  contact.  In  the  case  of  fric¬ 
tional  heat,  q  is  proportional  to  the  velocity  of  sliding,  v,  and  also  de¬ 
pends  on  the  total  force,  P,  exerted  between  the  two  relatively  sliding  ob¬ 
jects  and  the  number,  N,  of  contact  spots  at  which  they  touch.  For  plastic¬ 
ally  deformed  contact  spots  (which  is  by  far  the  most  frequent  case)  the 
following  expressions  are  found,  using  the  listed  symbols.  For  the  charac¬ 
teristic  velocity  vq  *  k^/ttNH^ / /?  =  p//P;  for  the  flash  temperature  in  gen¬ 
eral 

ATg  =  (M/PHg/TrN  v/A1)/(l/ZoS+Xr/So)  =  (pH^/c^) /(1/Z^+X^) 
where  W2>  =  l/(l+vr/3),  ZQ(vr>2)  =  (9/8)/(/T  +  1//8)  , 

S(vr<l)  -  Sq=  (4/3//(l+e3/4/3) (l+l/3e3/4)  and 

S(vr^l)  =  (e3^4+e3//4//8vr> /(l+e3^4/v',§vr) ; 

for  the  flash  temperature  at  very  low  speeds  for  circular  contact  spots 

AT  (v  «l)=u/PH  /nN  v/(X  +X  )  =  jUvP  v; 
or  s  i  i 

for  the  flash  temperature  at  very  high  speeds  for  circular  contact  spots 

ATo(vr>>1)=(9u/8Al)  (KlHsP/7TN)1/4  =  8|jpl/4 
and  for  the  flash  temperature  of  circular  contact  spots  at  v  =1 , 


(i) 


ATo(vr=l)  =  (uHg/p1c1)/(4/3+\r)  =  bp  . 

The  model  is  expanded  to  include  that  in  actual  fact  the  heat  is  evolved 
in  a  narrow  zone  in  the  softer  material  just  below  the  contact  spot.  For 
that  case  a  maximum  attainable  temperature  (if  below  the  melting  point)  is 
found,  of  roughly 

AT  =  p(10tt/3)H  /pc  =  fp>  attained  at  a  speed  in  the  order 
max  s  s  s 

of  lOOv  . 
o 

In  order  to  assist  in  the  practical  application  of  the  above  equations, 
tables  are  presented  of  approximate  values  of  the  parameters  p,  j,  g,  b  and 
f  for  a  number  of  metals  and  three  non-metals.  Pin/disk  asymmetries  arise 
because  at  high  speeds  the  heat  flow  is  practically  independent  of  the  pin 
material  but  is  determined  by  the  substrate,  unless  it  is  a  very  poor  heat  con¬ 
ductor  compared  to  the  pin. 

Lastly,  in  case  that  the  contact  spot  moves  relative  to  both  sides,  with 
speeds  v^  and  v^ ,  i.e.  relative  velocities  v^  =  k^R  and  v  2“  <2/R  wliere 
R  =  k'f/tt ,  it  is,  in  general, 

ATe  =  (p/PHs/7TN  v/X1)/[l/Zo(vrl)S(vrl)  +  VVvr2)S(v  J  ]  * 


'A 

1 

I 

List  of  Symbols 

'  « 

a 

Radius  of  circular  contact  spots 

b 

Parameter  for  the  computation  of  ATQ(vr=  1) 

B 

Workhardening  coefficient 

s 

c 

Specific  heat 

*  - 

e 

Ratio  of  axes  in  the  direction  of  sliding  and  at  right  angles  thereto 

*• 

of  elliptical  contact  spots 

1 

!- 

f 

Parameter  for  the  computation  of  AT 

r  max 

F 

Area  of  contact  spot 

. 

g 

Parameter  for  the  computation  of  ATo(vr»l) 

H 

s 

Hardness  of  the  softer  side 

j 

Parameter  for  the  computation  of  ATo(v^<<1) 

** 

L 

Sliding  distance 

/ 

m 

Width  of  the  deforming  zone  on  the  soft  side  in  units  of  a_ 

9 

n 

Measure  of  the  thickness  of  the  deforming  zone  on  the  soft  side 

9 

N 

Number  of  contact  spots 

P 

Parameter  for  the  computation  of  vq 

P 

Average  normal  pressure  at  contact  spots  (equal  to  Hg  in  case  of 

* ..  - 

plastically  deformed  contact  spots) . 

r.'  • 

r* 

P 

Normal  force  between  two  sliding  solids 

’ 

PN 

Normal  force  at  single  contact  spot 

PT 

Tangential  force  at  single  contact  spot 

-*■ 

• 

q 

Heat  input  per  unit  area  and  unit  time  at  contact  spots 

R  = 

yT7v 

Characteristic  dimension  of  contact  spots 

- 

; 

s 

Shape  factor  for  elliptical  spots  in  relative  motion 

• 

s 

O 

Shape  factor  for  elliptical  contact  spots  at  rest 

£  t 

O 

> 

II 

Characteristic  time  interval 

(iii) 

i# 

►  • 

.  .  v* 

.  '•  ■  .  •  •  . .  A;  '  '  '.-i 

T^  Melting  temperature 

ATq  Flash  temperature  of  circular  contact  spot  a'.  v^-value  indicated 

or  else  for  v  =  v 


AT  . 

ch 

AT 

e 

AT, 

n 

AT  . 

<, 

AT 


max 


v 

v 

crit 


v 

o 


v  = 

V* 


V 

W. 

in 


Y  =  v 


Z 

o 

^o 

Z 

<.0 

z 

z 

o 


£z 


o 

=  pHg/p^c^  Characteristic  temperature 

Flash  temperature  of  elliptical  contact  spots 
Flash  temperature  at  high  speeds 
Flash  temperature  at  low  speeds 

Maximum  temperature  within  the  deformed  zone  of  the  soft  side 
Velocity 

Speed  at  which  in  the  high-speed  approximation  AT^^  would  be 
reached 

k/R  Characteristic  velocity 
v/vq  Relative  velocity 
Heated  volume 

Mechanical  energy  input  into  one  contact  spot 
Rate  of  energy  input  to  one  contact  spot 
^Z^  Function  used  by  Jaeger  (2)  to  express  the  velocity  dependence 
of  flash  temperatures  in  case  of  friction  heat. 

Function  used  to  compute  the  velocity  dependence  of  flash  temp¬ 
eratures  of  circular  contact  spots 
Zq  at  high  speeds 
Zq  at  low  speeds 

Coordinate  normal  to  the  interface 

Depth  below  contact  spot  on  the  softer  side  at  which  shear  is  at 
maximum 

Depth  to  which  the  heat  diffusion  front  has  penetrated  into  the 
substrate 


(iv) 


Subsurface  shear  strain  on  the  soft  side 


Subsurface  shear  strain  at  z 

o 

Heat  diffusivity 
Coefficient  of  friction 

Subsurface  shear  stress  on  the  softer  side 

Critical  resolved  shear  stress  of  the  softer  material  at  small 
shear  strains 

Subsurface  shear  stress  on  the  softer  side  at  depth  z 


Subscripts 

Referring  to  substrate 
Referring  to  asperity  side 
Referring  to  number  of  contact  spots 

Referring  to  circular  spots, except  for  vq  which  is  applicable 
also  to  elliptical  spots. 


(  v) 


De-Mvstifying  Flash  Temperatures  -  II  Application  to  Heating  Through  Friction 


at  Plastically  Deformed  Contact  Spots 
by 

D.  Kuhlmann-Wilsdorf 
Department  of  Materials  Science 
University  of  Virginia,  Charlottesville,  VA  22901 

Introduction 

In  part  I  the  basic  equations  for  the  flash  temperatures,  AT,  of  ellip¬ 
tical  contact  spots  were  derived.  Invariably  AT  was  found  to  be  proportional 
to  q,  the  rate  of  heat  input  per  unit  area  of  contact  spot.  When  heating  is 
due  to  friction,  q  is  proportional  to  the  velocity.  Specifically,  for  a 
contact  spot  of  characteristic  radius  R  and  ellipticity  e  =  A/B,  with  A 
and  B  the  axes  in,  and  at  right  angles  to,  the  direction  of  sliding,  respect¬ 
ively,  we  have  in  accordance  with  Fig.  1,  with  R  =  /AB,  and  P„  the  normal  force 

N 

2  — 

q  =  v^v/ttR  =  ypv  ,  (la) 

where  p  is  the  average  normal  pressure  over  the  contact  area. 

The  value  of  p  depends  on  conditions.  When  is  gradually  raised  from 
zero  pressure,  p  increases  until  it  reaches  the  Meyer  hardness,  H,  (see  ref.  1) . 
Within  a  fairly  narrow  margin  that  pressure  is  then  maintained  until  the 
contact  spots  have  spread  to  the  point  that  they  occupy  a  sizeable  fraction 
of  the  macroscopic  contacting  area  and  they  begin  to  interact.  In  the  typical 
sliding  situation  the  actual  contact  area  is  much  smaller  than  the  macro¬ 
scopic  area  of  contact  and,  conversely,  the  normal  forces  exerted  are  large 
enough  to  cause  local  plastic  deformation  at  the  contact  spots  of  the  softer 
side.  In  these  cases,  then, 

q  =  UHgv  (lb) 

with  u  the  coefficient  of  friction  and  v  the  sliding  speed.  The  otherwise 


2 


redundant  subscript  s  is  used  here  as  a  reminder  that  H  is  the  hardness  of 

s 

the  softer  side,  whether  it  be  the  asperity  or  substrate  side. 

In  the  following  we  shall  examine  the  flash  temperatures  arising  in  this 
condition  of  plastically  deformed  contact  spots,  with  particular  emphasis 
on  unlubricated  "adhesive"  wear  between  metals. 

The  Basic  Equations 

a)  Circular  Contact  Spots 

Fur  unlubricated  "adhesive"  wear  with  plastically  deformed  contact 
spots,  the  relative  motion  between  the  two  sides  is  accomodated  through  sub¬ 
surface  shear  in  the  softer  side  according  to  Fig.  2.  We  shall  defer  discus¬ 
sion  of  the  consequences  of  the  fact  that  the  heat  is  developed  within 
either  the  asperity  side  or  the  substrate  side,  and  shall  for  the  time  being 
revert  to  the  model  of  Fig.  1  in  which  the  heat  is  evolved  in  the  interface. 
We  may  then  write  for  the  flash  temperature  of  a  circular  contact  spot  of 
radius  a  =  ^F/tt  -  R  and  area  F,  in  accordance  with  eq.  1  of  Part  I, 

ATo  -  yH^v/CXj/Z^Xp,  (2) 

where  the  symbols  are  mostly  the  same  as  in  part  I.  Specifically,  X  is  the 
thermal  conductivity  and  Zq  is  a  function  of  relative  velocity,  v  , 

vr  =  v/vo  =  v/ (k^/R)  (3) 

where 

vo  =  (k^/R)  =  X1/p1c1R  .  (4) 

Here  is  the  thermal  dif f usivity, p  is  the  mechanical  density, and  c^  is 

the  specific  heat  of  the  substrate  side  ,  -  meaning  the  side  relative  to 
which  the  contact  spot  moves,  while  it  is  stationary  on  the  asperity  side. 

As  was  shown  in  part  I, 

W2)  :  Z!o  ■  1/(1+V3> 


(3a) 
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Z  (v  >2)  =  Zu  =  (9/8) /(.v  +l/v b) 
or  ho  r 


which  for  quite  large  values  of  v^_  reduces  to 


Z  (v  »1)  =  9/8v'v 
o  r  r 


Introducing  the  relative  thermal  conductivity  of  the  asperity  side  as 

Ar  =  X2/\,  (f 

the  characteristic  temperature  rise  as 

“Tch  ■  “HsRv<A  ■  pHs/plV  <; 

and  the  function 

Y  =  v  Z  ,  (6 

r  o 

ao  used  by  Jaeger  (2),  we  can  rewrite  eq .  2  into  the  more  convenient  form 


AT  =  AT  /(1/Y+A  /v  )  =  AT  v  /(1/Z+a  ). 
o  ch  rr  chr  or 


Fig.  3  shows  AT  /AT  ,  as  a  function  of  A  and  v  ,  using  Z  according 
6  och  rr°o  ° 

to  Jaeger's  computation  (Fig.  7  of  ref.  2)  as  reproduced  in  Fig.  2  of  Part  I. 

As  seen,  for  small  speeds,  when  Z  =1,AT  /AT  ,  rises  linearly,  namelv  as 

o  o  ch 

(AT  /AT  )  =  v  /(1+A  )  (9b) 

o  ch  v<<v  r  r 
o 

and  at  intermediate  and  high  speeds  when  eq.  5b  applies, 

(AT  /AT  .  )  =  v  /[ (8/9) (/v  +l//8)+A  ]  (9c) 

o  ch  vA2v  r  r  r 

o 

in  accordance  with  the  corresponding  expressions  for  Zq  (eqs.  5b  and  5c). 
b)  Elliptical  Contact  Spots 

ror  elliptical  contact  spots  the  appropriate  expression  for  the  flash 
temperature  is,  corresponding  to  eq.  35  of  Part  I, 

1Te/1Tch  ■  im/YS+ySovt)  -  vr/(l/ZoS+>r/So>  . 


The  shape  factor  S  was  derived  in  Part  I.  For  vr^l  it  is 


(10) 


(11) 


S(v  >1)  =  (e^Ve^^/vSv  ) / ( 1+e 3 ^ /v 8v  ) 
r  r  r 

and  for  v  <1  it  is 
r 

S(vr<l)  s  Sq  =  (4/3)/(l+e3/4/3)(l+l/3e3/4)  .  (12] 

Correspondingly,  for  velocities  below  v  ,  Sq  appears  in  eq.  10  not  only  in 

coniunction  with  A  but  also  with  Z  so  that 
r  o 

(AT  /IT  )  =  S  v  / (1/Z  +\  )  .  (13; 

e  ch  v  11  or  or 

r 

In  fact,  except  for  the  introduction  of  AT  ,  =  iiH  /d,  c,  in  lieu  of 

ch  s  1  1 

q v  F / rr' /  ,  and  multiplication  with  v  this  is  eq.  3  4  of  Part  I.  Figs.  4  to  6 

show  IT  /AT  ,  for  the  same  choice  of  parameters  for  which  AT  /(qvF/ir/X  ) 
e  ch  o  1 

was  computed  and  represented  in  Figs  7-9  of  Part  1. 


Application  to  “Adhesive"  Sliding 


a)  The  Model 


While  Fig.  1  is  an  oversimplification  of  the  actual,  vastly  more  compli¬ 
cated  case  of  frictional  heating  at  contact  spots  between  sliding  solids, 

Fig.  2  is  still  an  idealization,  although  rather  more  realistic.  It  envisage, 
that  the  mechanical  work  input 

W.  =  yP  (14) 

in 

when  two  solids  slide  past  each  other  for  unit  distance  under  the  normal 
force  P  is  totally  converted  into  heat  via  plastic  deformation.  Indeed, 


the  fraction  of  which  is  stored  in  the  form  of  crystal  defects  is  typ¬ 

ically  very  small  (e.g.  ref.  3)  and  it  may  well  be  overbalanced  by  additional 
heating  through  chemical  reactions  at  the  surface  such  as  through  oxidation. 

In  agreement  with  a  previously  developed  model  (4,5),  it  is  assumed 
that  the  shearing  takes  place  in  the  softer  material  underneath  the  contact 
spots  and  ls  essentially  confined  to  a  zone  ma  wide  and  a/n  thick  where 


5 


2<m<4  and  2Tn<50  with  m~3  and  n'10  good  average  values. 

Since  according  Co  this  model  the  shearing  zone  typically  is  much  thicker 
than  incidental  surface  films,  their  presence  will  not  greatly  affect  the 
results.  However,  it  has  already  been  noted  that,  strictly  speaking,  the 
heat  evolution  does  not  take  place  at  the  interface,  where  indeed  there  will 
normally  be  relative  rest  and  thus  no  heat  evolution  of  any  significance  at 
all,  but  within  the  softer  side  quite  close  to  the  interface. 

This  concept  is  akin  to  the  proposal  of  Rozeanu  and  Pnueli  (6,7)  that  shear 
instabilities  could  concentrate  the  relative  displacement  into  a  narrow 
subsurface  region  which  becomes  strongly  heated.  Yet,  the  present  model  has 
a  different  physical  basis,  since  there  is  always  enforced  relative  rest 
at  the  interface  on  account  of  micro  interlocking  to  a  relief  depth  of,  say, 


z^  on  average,  where  zo$1000A(ref .  4).  For  the  remainder  at  depths 

z >z  the  shear  stress  appears  (4,8),  to  drop  as 
T  =  TsZo/‘,(z+Zo), 


(15) 


lie  a  representative  work  hardening  law,  linking  subsurface  shear,  7 ,  to 
the  prevailing  stress,  T,  is  (4) 


=  T  +  B  ; . 


(16) 


Consequently,  after  repeated  passes, or  a  sufficiently  long  sliding  distance 
that  quasi-equilibrium  is  established,  we  expect 


Y  =  Y  exp-  :  (z  -z)/B(z  +z) ;  (17) 

s  s  o  o 

wi  th 

Y  =  exp ■('-') /B;  (18) 

s  so 

Previously,  (  -»  )  values  for  T  ,  :  and  B  were  semi-quantitativelv  in- 

o  s 

f erred  fr-tn  published  measurements  of  -  for  a  worn  copper  pin  (9  )  and  were 
found  to  he  compatible  with  the  known  materials  properties  and  sliding 
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MATERIAL 

p/]03 

c 

A 

* /IQ'4 

H/108 

p19 

B* 

G10 

F* 

th 

Ag 

10.5 

233 

418 

1.71 

3-7 

21 

88 

50 

2100 

960°C 

Al 

2.7 

380 

210 

0.86 

1.4-4 

8.2 

52 

47 

1300 

660 

Au 

19.3 

129 

310 

1.18 

2-7 

14.0 

77 

51 

1900 

1063 

Be 

1.85 

1750 

150 

0.46 

15 

10.0 

200 

160 

4900 

1280 

Cr 

7.1 

440 

67 

2.17 

7-13 

3.8 

140 

186 

3400 

1903 

Ci, 

8.9 

381 

380 

1.14 

4-7 

15 

70 

48 

1700 

1083 

Fe 

7.8 

452 

60 

0.17 

6 

2.3 

73 

125 

1300 

1539 

Ir 

22.4 

134 

50 

0.167 

27 

4.9 

380 

460 

9400 

2443 

Mg 

1.74 

1005 

159 

0.91 

2 

7.2 

49 

48 

1200 

657 

Mo 

10.0 

250 

140 

0.56 

18 

13 

310 

222 

7500 

2610 

Ni 

8.8 

446 

70 

0.178 

7-20 

3.7 

150 

209 

3600 

1453 

Pb 

11.3 

126 

35 

0.25 

0.5 

1.0 

15 

40 

370 

327 

Pd 

12.2 

247 

70 

0.25 

4-10 

3.6 

100 

143 

2500 

1552 

Pt 

21.4 

132 

70 

0.25 

4  -  8 

3.4 

91 

130 

2200 

1772 

Sn 

7.3 

223 

64 

0.38 

.45-. 6 

1.5 

14 

28 

337 

232 

Ta 

16.6 

138 

54 

0.24 

12 

4.7 

220 

280 

5500 

2996 

Ti 

4.5 

611 

17 

0.062 

11 

1.2 

170 

420 

4200 

1668 

W 

19.3 

134 

180 

0.70 

12  -  40  20 

430 

250 

10,500 

3380 

Zn 

7.1 

384 

111 

0.41 

3  -  6 

4.9 

71 

85 

1700 

420 

Brass 

8.5 

377 

50 

0.33 

9 

5.6 

120 

230 

2900 

Steel 

7.8 

490 

50 

0.12 

15-30 

3.2 

250 

340 

6200 

Glass 

2.4 

840 

1.2 

.0058 

8 

0.09 

170 

1400 

4200 

Brick 

2.6 

840 

0.8 

.0038 

40  (?) 

0.13 

780 

5800 

19,000 

Gran ite 

2.4 

950 

1.7 

0.012 

30 

0.37 

560 

3900 

14,000 

TABLE  I 
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by  a  factor  of  3  for  Ni/Ti,  of  ^10  for  steel/granite,  and  ^6  for  Mg/Cu. 
The  values  in  Table  III  suggest  that  flintlock  ignition  is  theoretically 
the  best,  indeed  the  only  realistically  feasible,  manner  to  generate  a 
spark  in  one  strike. 

13)  Pb  and  Sn  are  much  softer  than  engineering  materials.  Hence  these  deform 
when  paired  with  those,  and  the  flash  temperatures  are  correspondingly 
low.  Thus  Pb  and  Sn  platings  provide  fine  lubrication  layers. 

14)  Hard  platings  by  contrast  can  transfer  the  shear  strain  to  the  other 
side  that  may  be  less  sensitive  or  less  used,  e.g.  from  pin  to  disk. 

Also,  hard  platings  will  mostly  reduce  p  and  thus  the  flash  temperatures. 

15)  The  theory  may  be  readily  expanded  to  the  case  in  which  the  contact 
spot  moves  relative  to  both  sides,  as  may  occur  in  gears.  In  that  case. 


in  terms  of  the  contact  velocities  v  and  v 2  relative  to  the  two  sides, 
with  v  =  v^  +  v^,  it  is  v^  =  k^/R,  v^  =  <^/R  and 


-T  / (i* PH  j 
e  s 


™  *>  ’  1/1  Wvrl>S<Vrl) 


^Zo^Vr2^Vr2^  ^ 
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coefficient  of  friction,  -  in  actual  fact  thus  on  the  shear  strength 
of  the  softer  material. 

7)  The  data  in  Table  1  yield  this  maximum  temperature  to  be  above  the 
melting  point  if  the  coefficient  of  friction  1  '  es  above  "■  1/2,  on 
average.  However,  in  fact  the  melting  temperature  of  the  softer  side 
cannot  be  much  exceeded  because  the  material  loses  its  shear  strength 
and  thus  the  coefficient  of  friction  drops  strongly 

8)  The  critical  velocity  at  which  the  maximum  temperature  is  reached  is 

estimated  at  about  100  v  or  moderately  less. 

o 

9)  Again,  using  characteristic  values  of  the  parameters  involved,  vq  is 
found  as  low  as  'VLm/sec  at  P=1N  for  lead,  and  as  high  as  n^Om/sec  for 
silver  and  for  tungsten.  Since  vq  is  proportional  to  1//P,  it  can  be 
as  low  as  lem/sec  for  railroad  wheels. 

10)  The  flash  temperature  at  v=vq  is  independent  of  P,  but  is  proportional 
to  u •  For  y-1/3  it  ranges  between  ^5°C  for  Sn  and  140°C  for  W. 

11)  For  low  speeds,  ATq  = j  u/p  v  where  j depends  on  both  of  the  materials. 
For  Sn/Pb,  for  example,  j~13  and  for  Ag/Sn  it  is  j~3  in  mks  units. 
However,  j  is  quite  high  for  hard  materials  and  those  with  low  thermal 
dif fusivities ,  (e.g.  it  is  js90  for  titanium  on  iridium). 

12)  In  the  high-speed  regime,  pin/disk  asymmetries  arise.  These  can  be 
further  accentuated  because  geometry  and  heat  flow  favor  long-lived  de¬ 
forming  asperities  on  the  pin  and  therefore  the  pin  tends  to  assume 

the  role  of  the  soft  material  even  if  ordinarily  it  would  be  the  harder 
one.  Therefore,  while  the  disk  controls  the  heat  flow, the  pin  typically 
controls  the  heat  generation.  A  table  of  estimated  values  shows  that 
on  reversal  of  pin/disk  roles  for  the  same  materials,  the  temperature 
rise  at  the  interface  can  differ  by  a  factor  of  several;  specificallv 
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lubricancs.  Conversely,  Ti  tends  to  run  "hot"  among  the  metals,  both  on 
account  at  high  hardness  and  of  low  heat  conductivity. 

The  by  far  lower  heat  conductivities  of  ceramics  as  compared  to  metals, 
finally,  cause  flash  temperatures  of  these  to  be  very  much  higher,  to  which 
their  generally  greater  hardness  still  adds.  Interestingly,  according  to 
table  III,  trial  and  error  led  to  the  best  choice  for  making  fire  in  striking 
flint  against  steel  in  flintlock  rifles. 

Summary 

1)  Flash  temperatures  cannot  be  given  accurately  because,  in  nature,  they 
are  ill-defined  and  vary  greatly  from  place  to  place  and  moment  to  moment. 

2)  Even  so,  good  estimates  are  possible,  provided  essential  parameters 
are  known. 

3)  Speeds  are  slow  or  fast  relative  to  a  characteristic  speed  vq.  At  that 
speed  the  "diffusion  distance"  of  the  heat  energy  equals  the  character¬ 
istic  dimension  of  the  contact,  a,  within  the  time  6t  during  which  the 
substrate  moves  through  the  length  of  its  radius  -  more  or  less.  I.e. 

in  the  interval  dt  s  a/v  the  diffusion  distance  is,  Sz  =  /<6t  =  a.  while 

o  * 

the  spot  moves  through  a  =  /<a/vQ  . 

4)  Significantly  below  v  ,  the  asperity  acts  much  as  if  it  were  at  rest, 
and  the  heat  flows  into  the  two  materials  in  proportion  to  their  heat 
conductivities,  X^  and  X^ . 

2 

5)  At  high  speeds,  and  provided  that  vr>>^r  >  the  softer  material  heats 
adiabatically  through  rapid  shear  at  the  interface,  but  the  heat  flows 
practically  exclusively  into  the  substrate,  so  that  only  its  heat  con¬ 
ductivity  enters  into  the  problem. 

6)  In  the  fully  adiabatic  case,  when  shearing  takes  place  only  in  the  softer 
material,  a  maximum  temperature  is  reached  which  depends  on  the  size  of 
the  shear  zone  relative  to  the  size  of  the  contact  spot  and  on  the 


answers  presented  in  this  paper  are  upper  limits  as  far  as  the  discussed 
shift  of  the  zone  of  heat  evolution  is  concerned. 

b)  Reliability  of  Results 

It  will  have  become  clear  that  the  critical  parameter  in  calculations 
of  flash  temperatures  is  vq  since  it  is  in  comparison  to  this  velocity  that 
actual  velocities  must  be  judged  as  "slow"  or  "fast".  Fortunately,  vq  de¬ 
pends  on  the  root  of  both  the  number  of  contact  spots  and  hardness,  so  that 
errors  in  these  parameters  do  not  affect  the  result  sensitively.  However, 
as  seen  in  table  I,  for  same  P  and  N,  the  characteristic  velocity  v°  has  very 
widely  different  values  depending  on  the  material,  while  it  decreases  with 
increasing  P,  and  increases  when  N  increases. 

If  the  input  parameters  are  all  known,eq.  10  in  conjunction  with  eqs . 

5a  and  5b  for  the  determination  of  Z  and  eqs.  11  and  12  for  S  and  S  ,  re- 
spectively,  will  give  results  within  a  few  percent  of  the  accurate  mathe¬ 
matical  solution  for  model  I.  This  is  amply  satisfactory  and  no  further 
refinements  appear  to  be  worthwhile  on  account  of  the  unavoidably  larger  errors 
and  uncertainties  which  arise  due  to  discrepancies  between  the  model  and 
actuality,  e.g.  as  in  Fig.  2  and  through  erratic  changes  in  N  and  spot  shape. 

Further,  a  more  thorough  study  will  have  to  include  the  changes  of  X,  p,  < , 

c  and  H  with  temperature.  Fortunately,  with  the  exception  of  H  ,  these  are 
s  s 

minor . 

c)  Qualitative  Features  of  Results 

Notwithstanding  the  above  limitations,  the  results  obtained  are  very 
helpful.  Clearly,  for  example,  under  otherwise  same  conditions  materials  run 
"hot"  or  "cool"  principally  depending  on  their  thermal  conductivity  and 
hardness.  Thus,  as  evident  from  Tables  I  to  III,  Pb  and  Sn  always  run  "cool" 

on  account  of  their  very  low  hardness,  and  thus  they  may  be  used  as  effective 


Lastly  Table  III  demonstrates  the  pin/disk  asymmetry  for  a  number  of 
material  pairings,  assuming  for  the  pin  either  one  half  of  the  listed  hard¬ 
ness,  _or  the  lowest  value  in  the  range  of  listed  for  the  pin  material. 

No  values  of  j  are  given  if  those  pin  hardnesses  lie  above  the  upper  hard¬ 
ness  value  given  for  the  substrate,  on  the  assumption  that  in  such  a  case 
the  substrate  material  will  continue  to  slide.  As  already  indicated,  pin/disk 
asymmetry  will  be  also  expected  in  that  case,  but  the  corresponding  values 
can  be  readily  calculated  and  might  confuse  the  reader  if  they  were  also  in¬ 
cluded  in  Table  III. 

Discussion 

a)  Effect  of  Sub-Surface  Instead  of  Interfacial  Heat  Evolution 

The  above  graphs  (Figs.  3  to  6)  and  Tables  I  to  III  should  go  far  to 
aid  experimentalists  in  making  rapid  estimates  and/or  calculations  of  the  flash 
temperatures  due  to  friction  in  unlubricated  adhesive  sliding.  With  some 
caution  these  may  be  also  used  in  the  presence  of  minor  surface  films,  but 
not  necessarily  if  sliding  is  not  "adhesive".  The  deciding  factor  is  whether 
the  heat  is  evolved  at  restricted  contact  spots  sufficiently  small  and  suf¬ 
ficiently  widely  separated  that  the  models  of  Fig.  1  or  2  are  applicable. 

As  to  the  displacement  of  the  zone  of  heat  evolution  somewhat  into  the 
interior  of  the  softer  material,  this  will  have  the  lesser  effect,  the  more 
closely  to  the  surface  it  is  situated.  At  any  rate,  the  highest  temperatures 
will  be  in  the  softer  material  and  the  heat  flow  into  the  harder  side  will 
be  diminished.  Correspondingly,  also,  the  surface  itself  will  be  at  a  lower 
temperature  than  calculated  on  the  basis  of  the  model  in  Fig.  1. 

By  and  large  it  appears,  intuitively,  that  the  above  effects  will  not  be 
major,  and  for  this  reason  the  theory  has  been  presented  making  use  of  the  model 
in  Fig.  1,  and  the  contact  spot  size  has  been  used  rather  than  the  somewhat  larger 
area  of  actual  heat  evolution.  Future  investigations  will  have  to  determine 


the  magnitude  of  the  consequent  error.  At  any  rate,  it  would  seem  that  the 
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Further,  in  computing  p^,  b*,  g  and  f*  in  table  I,  the  average 
value  of  H  was  used  where  a  range  rather  than  a  single  value  was  given. 

At  any  rate,  hardnesses  vary  widely  from  specimen  to  specimen,  as  also 
with  temperature,  and  therefore  all  of  the  values  given  for  the  computed 
parameters  are  guide-line  values  only,  and  certainly  not  accurate.  The  use 
of  those  parameters  is  in  accordance  with  eqs.  22,  25,  26  and  29  namely, 
for  circular  spots, 

(i)  Find  v  ,  the  characteristic  velocity  for  the  listed  material  on  the 
substrate  side,  as 

Vo  =  P10//?  (33) 

(ii)  Find  the  flash  temperature  for  the  material  sliding  on  itself  with 

speed  v  as 
o 

ATQ(vr=l)  =  ub*  (34) 

(iii)  Find  the  flash  temperature  in  the  case  of  /v^>>Ar,  and  the  listed 
material  as  the  substrate,  as 

ATo(vr>>A^)  =  g10l//p  Fv  (35) 

(iv)  Find  the  estimated  maximum  temperature  in  the  listed  material  when  it 
acts  as  the  softer  side  as 

AT  =  fy  (36) 

max 

provided  that  ^max<^ »  since  otherwise  melting  inhibits  further  rise  of 
energy  input.  Note,  therefore,  that  the  ratio  ^T  ^/f  provides  an  estimate 
of  the  lowest  coefficient  of  friction  which  can  lead  to  melting. 

In  table  II  are  listed  values  of  j  according  to  eq.  24.  With  their 
aid  the  flash  temperatures  at  low  speed  (i.e.  v  <1)  may  be  found  as 

AT  (vr<l)  =  y/p  v  j  .  (37) 

Since  low  speeds  are  concerned  here,  there  is  no  pin/disk  asymmetry.  Again 
the  average  values  of  are  chosen,  extracted  from  Table  I. 


exists  for  the  simple  reason  that  on  reversal  of  the  two  materials  A  changes 
from,  say,  1/2  to  2,  or  4  to  1/4,  etc.  As  seen  from  Fig.  3  this  causes  the 
corresponding,  typically  quite  marked  change  of  flash  temperatures.  To  the 
extent  that  wear  rates  go  up  with  temperature  this  means  that  making  the 
pin  of  the  more  highly  heat  conducting  material  will  generally  decrease  the 
wear  rate. 

It  may  be  added  here  that  the  time  intervals  in  which  contact  spots  are 
stationary  on  the  substrate  side  need  be  no  longer  than  2a /v  or  /eF / tt'/v , 
respectively,  to  make  the  equations  applicable.  If  the  spot  moves  in  regard 
to  both  sides,  the  appropriate  values  of  Zq  and  S  must  be  applied  also  in 
connection  with  A^  in  all  of  the  above  equations. 

Numerical  Examples 

The  experimentalist  is  often  interested  in  some  rough  estimates  but  not 
interested  to  go  through  the  equations  in  detail.  In  such  cases  the  curves 
in  Figs.  3-6  should  be  very  helpful.  To  simplify  the  task  even  more,  Tables 
I  to  III  present  some  approximate  values  for  the  hardness  and  the  parameters 
p  (eq.  22b),  j  (eq.  24b),  b  (eq.  25b),  g  (eq.  26b)  and  f  (eq.  29b). 

Specifically,  in  Table  I  are  collected  the  values  for  p,  c.  A,  K,  H,  p, 
b,  g,  f  and  the  melting  temperatures  for  a  number  of  materials,  -  in  regard 
to  specific  conditions  which  are  explained  in  the  legend  to  Table  I,  where 
also  the  significance  and  use  of  those  parameters  is  more  fully  explained. 

Throughout,  N  =  10,  m  =  3,  and  n  =  10  have  been  assumed.  It  is,  of 
course,  a  simple  task  to  convert  to  other  values  of  N,  m  and  n;  e.g.  indi¬ 
cating  the  chosen  value  of  N,  the  number  of  contact  spots,  by  the  corresponding 
subscript  (as  done  in  Table  I) ,  we  have 

PN  =  P1()vWl0,  (32) 

and  similarly  in  regard  to  m  and  n  whose  reasonable  ranges  are  believed 


to  be  2  <  m  £  4  and  4  <  n  <  50  while  1  <  N  $  100. 
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Pin-Disk  Asymmetries 

The  above  expressions  for  vcri(./vo  and  ATmax/AT0 (vr=l)  are  somewhat 

cumbersome  on  account  of  the  factor  p, c, /p  c  .  This  reduces  to  unitv  if 

1  1  s  s  J 

the  substrate  is  the  softer  material  but  otherwise  can  lie  within  a  rather 
wide  range  of  values  depending  on  the  materials  concerned.  At  any  rate, 
Vcr-Lt/V0  wiH  almost  certainly  be  much  larger  than  unity  so  that  the  high 
speed  approximation  for  Zq  applies,  unless  is  rather  large  (see  Fig.  3) 
in  which  case  eq.  2  or  10  in  conjunction  with  eq.  5b  must  be  applied. 

At  first  glance  it  may  appear  strange  that  the  "substrate"  side  might 
be  the  softer  side  just  as  well  as  the  "asperity"  side,  for  the  reason  that 
one  intuitively  tends  to  think  of  the  asperity  as  being  sheared.  However, 
on  which  side  the  contact  is  stationary,  whether  for  short  or  extended  periods, 
is  not  so  much  a  question  of  materials  properties  than  of  macroscopic  geometry: 
If  one  of  the  two  members  is  much  smaller  than  the  other,  such  as  in  a  pin 
and  disk  or  block  and  ring  arrangement,  contact  spots  can  be  stationary  only 
on  the  smaller  object,  i.e.  the  pin  or  disk,  but  not  on  the  larger  one. 

A  stationary  contact  spot  tends  to  bulge  out  through  thermal  expansion, 
which  increases  the  load.  This  then  causes  further  expansion  in  the  typical 
thermo-mechanical  behavior  (compare  ref.  10).  Further,  the  bulk  temperature 
of  the  smaller  member,  i.e.  normally  the  asperity  side,  will  tend  to  rise  more 
and,  lastly,  as  already  discussed,  the  heat  is  not  strictly  speaking  generated 
in  the  interface  but  on  the  softer  side.  All  of  these  effects  have  the 
tendency  to  favor  sliding  in  the  smaller  member,  i.e.  the  pin,  and  thus  cause 
its  wear,  even  if  normally  its  material  should  be  the  harder. 

This,  then,  is  a  cause  of  strong  pin/disk  asymmetry.  However,  even  if 
the  discussed  reversal  in  the  role  of  the  "softer"  material  does  not  occur 
on  reversal  of  pin  and  disk  materials,  yet  a  significant  pin/disk  asymmetry 
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the  width  of  which  is  controlled  geometrically.  Let  us  then  ask  what  the 

maximum  temperature  is  which  may  be  attained  in  that  zone.  To  this  end 

we  proceed  much  as  in  deriving  the  equations  in  Part  I,  as  follows: 

The  enert y  input  when  a  contact  spot  moves  through  distance  L  is 

W.  =  y(P/N)L  =  p7ra2H  L  (27) 

in  s 

which  is  expended  in  the  volume 

v  =  ma(a/n)L  (28) 

according  to  Fig.  2.  Provided  that  the  velocity  is  high  enough  so  that  the 
heating  in  that  volume  is  essentially  adiabatic,  the  resulting  temperature 
rise  is 


AT  =  W.  /p  c  V  =  yirH  n/p  c  m  =  fy 
max  mss  s  s  s 


(29a) 


f  =  ttH  n/p  c  m 
s  s  s 


(29b) 


where  the  subscript  s  reminds  us  that  the  softer  material  is  concerned. 

Comparing  eq.  29  with  eqs.  22,  25  and  26  shows  that  the  above  maximum 
temperature  would  be  reached  in  the  high-speed  case  at  the  critical  speed 

Vcrit  given 

v  _,../v  =  [  (87Tn/9m)  (p  c  /p  c  )]2  (30) 

crit  o  i  1  s  s 

and  that  AT  relates  to  AT(v  =  1)  as 
max  r 

AT  /AT  (v  =  1)  =(TTn/m)  (4/3+X  /A  )  (p.  c  /p  c  )  .  (31) 

max  or  2  1  1  1  s  s 

Physically,  the  significance  of  AT  is  this:  Provided  that  the  values 

max 

of  AT  computed  according  to  eqs.  23  or  26a  lie  well  below  AT^^,  they  are 
presumed  to  be  rather  reliable  and  the  heatflow  to  be  well  described  by  the 


model  of  Fig.  1.  However,  as  the  computed  values  approach  AT  ,  they  have 

max 

to  be  regarded  with  skepticism  and  are  in  all  likelihood  too  large. 


(22a) 


with 


v  =  K.  /ttNH  /7p=p//P 
o  1  s 


p  =  k1v/ttNHs  (22b) 

and  the  flash  temperatures  for  circular  contact  spots  (eqs.  2  and  9)  are 
found  as 

ATq  =  y/PH  /irNv/X  (1/z  +X  )  (23) 

in  general.  At  very  low  speeds  eq.  23  reduces  to 

AT  -  u^»/(X  +X  )  =  jy/p  v  (24a) 

with 

j  =  /Hs7ttN  /(X1+X2)  ,  (24b) 

while  at  the  characteristic  speed  of  v  it  is 

o 

ATQ(vr=l)  =  (yHs/plCl)/(4/3  +  X  )  =  by  (25a) 

with 

b  =  (Hs/p1c1)/(4/3+Ar)>  (25b) 

and  at  high  speeds  such  that  k\T>>X  and  vr>2 ,  at  the  least, 

ATQ(vr»l)  =  ATh  =  (9y/8A1)v^(K^P/TrN)1/4  =  gyP1/4  &  (26a) 

with 

g  =  (9/8a1)(k^/ttN)1/4  .  (26b) 


The  Maximum  Attainable  Temperature 

Eq.  26  can  clearly  not  apply  to  any  arbitrarily  high  velocity.  The  most 
obvious  limit  to  the  flash  temperature  is  the  melting  point  of  either  of 
the  two  sides.  Once  this  is  attained,  the  coefficient  of  friction  must  drop 
preciptously  (ice-skating  being  the  best-known  example)  and  thereby  the 
work-input . 

According  to  Rozeanu  and  Pnuely  (6,  7)  melting  can  occur  through  a 
thermo-mechanical  instability.  This  may  well  be  possible.  However,  as  was 
discussed  above  in  connection  with  Fig.  2,  one  must  expect  the  heat  genera¬ 
tion  not  to  occur  at  the  interface  in  any  event,  but  within  a  narrow  zone 
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conditions.  The  expended  mechanical  energy  density  found  from  15  to  18  is, 


for  z>z 


x  Y  =  T  y  [ 2z  /(z  +z) ] exp{ (x  /B) (z  -z)/(z  +z)} 
s  s  o  o  s  o  o 


Since  invariably  (x  / B)  will  be  larger  than  unity  and  more  typically  be  in 
s 

the  order  of  10,  the  decrease  of  energy  density  input  with  depth  described 

by  eq.  19  is  very  rapid.  E.g.  for  x^/B  =  10  it  decreases  by  a  factor  of  40 

between  z  =  z  and  z  =  2z  . 

o  o 

Thus  the  zone  of  heat  evolution  will  be  very  narrow  and  hence  be  well 
approximated  by  a  mathematical  plane  as  in  Fig.  1  for  the  great  majority  of 
cases.  That,  even  so,  the  parameter  n  appearing  in  Fig.  2  is  not  estimated 
to  be  larger  than  50  and  more  typically  is  believed  to  be  about  10,  is  for 
reasons  of  geometry,  in  that  contact  spots  cannot  be  expected  to  be  accurately 
Planar  over  their  whole  extent  but  will  be  warped.  Indeed,  one  may  feel 


rather  comfortable  with  the  estimate  of  n  =  10,  being  confident  that  it  will 
mostly  be  good  to  within  a  factor  of  less  than  two. 
b)  The  Parameters  q  and  a  (or  R) 

Overwhelmingly,  the  softer  side  is  plastically  deformed  and  its  com¬ 
pression  hardness,  Hg,  determines  the  total  contacting  area  under  total 
normal  load  P.  Therefore,  if  there  are  N  contact  spots,  assumed  to  be  cir¬ 
cular,  their  radius  is  found  as 

a  =  /P/ttNHs  .  (20) 

2 

If  now,  ->ver  the  total  contact  area  NTTa  ,  the  frictional  work  is  dissipated  at 
the  rate  W_^  =  pPv,  the  power  density  q  becomes 


q  =  W.  /Nxa  =  uH  v 
in  s 


(see  eq .  lb ) . 


With  this,  the  characteristic  velocity  (eq.  4)  is  found  as 


i 
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eqs .  10  co  12  for  X  =0. 

r 

Fig.  5  As  Fig.  4  but  for  A  =1. 
Fig .  6  As  Fig.  4  but  for  X^=10. 
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Table  L: 

Approximate  values,  in  mks  units  throughout,  of  density  (p)  ,  specific 
heat  (c) ,  thermal  conductivity  (A),  thermal  diffusivity  (k)  ,  hardness  (H) , 
and  melting  temperature  (T  )  for  a  number  of  materials,  together  with  the 
parameters  p^  (i.e.  p  of  eq.  22b  for  the  case  of  N  =  10  and  the  material 
serving  as  substrate),  b*  (i.e.  b  of  eq.  25b  for  the  case  of  the  two  sides 
being  of  the  same  material),  g  (i.e.  g  of  eq.  26b  for  the  case  of  N  =  10, 
and  the  material  serving  as  substrate  and  being  the  softer  member) ,  and  f* 

(i.e.  f  of  eq.  29b  for  the  case  that  the  listed  material  serves  as  the  softer 
member,  that  n  =  10,  and  that  m  =  3) . 

The  values  of  those  parameters  and  their  use  (for  circular  contact  spots 
under  the  named  conditions  and  choosing  the  average  value  of  H  where  a  range  is 
given)  are  as  follows: 

P^q  =  <»10nH  ;  vq  -  p/^P 

b*  =  (H/pc) / ( 4 / 3  +  1)  ;  AT  (v  -1)  =  ub 

o  r 

g10  =  (9/8A)(k2H3/10-)1/4  ;  ATo(vr»l)  =  gyP1/4v^ 

f*  =  10ttH/3pc  ;  AT  =  fy 

max 

Thus  p^Q  gives  an  indication  of  the  characteristic  speed.  E.g.  under 
otherwise  same  conditions,  is  roughly  20  times  larger  for  silver  than  for  lead, 
and  200  times  larger  than  for  glass. 

Further,  gives  a  qualitative  impression  whether  the  respective  materials 
tend  to  "run"  hot  or  cool  at  high  relative  speeds.  Namely,  large  values 
indicate  relatively  high  flash  temperatures,  i.e.  hot  contact  spots, and  low  values 
cool  contact  spots.  Hence,  tin  runs  coolest  (and  therefore  is  a  good  lubricant 
plating)  while  the  non-metals  run  hot. 

Next  f*  is  useful  in  judging  (roughly)  whether  a  material  is  liable  to  melt, 
i:  r  undergo  some  phase  change  etc.,  through  friction,  since  y  -  ^  x/f  is  the 
lowest  coefficient  of  friction  at  which  some  specific  temperature  AT^^  (e.g. 
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that  at  which  austenite  transforms  to  martensite)  can  be  reached. 

Lastly,  b*  is  a  guide  to  the  flash  temperature  that  is  to  be  expected 
at  intermediate  speeds  for  any  given  value  of  y. 

Note  that  AT  is  proportional  to  the  coefficient  of  friction  in  all  cases. 
Hence  the  most  obvious  corrective  measure  in  order  to  decrease  flash  temp¬ 
eratures  would  be  to  decrease  y,  e.g.  through  lubrication. 
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Table  II 

Values  of  =  vH^/ttN/CX^  +  X  )  according  to  eq.  24,  for  N  =  10, 
based  on  the  values  in  Table  I,  assuming  the  average  value  of  H  where  a 
range  is  given.  From  these  find  the  temperature  rise  for  circular  contact 
spots  at  low  speeds;  -  for  v^  <  0.7  with  an  error  of  less  than  20%,  for 
v^_  <  3  within  better  than  a  factor  of  two,  and  in  each  case  obtaining  an 
upper  limit.  Namely,  AT^  =  j v.  Since  under  otherwise  same  sliding 
conditions  the  flash  temperature  is  proportional  to  j ,  as  seen,  it  is 
evident  how  much  hotter  the  contact  spots  are  for  non-metal  on  non-metal  than 
when  at  least  one  metal  is  involved.  Note  the  consistently  low  values  of  j 
when  at  least  one  of  the  sliding  member  is  a  very  good  heat  conductor,  e.g. 
Ag,  A1 ,  Au,  Cu  and  Mg,  while  values  involving  Sn  are  low  principally  on 
account  of  its  low  hardness. 


25 


Table  III 

Selected  g^values  as  in  Table  I  with  N  =  10,  but  (i)  choosing  Hg  as  either 
one  half  of  the  H  value  given  in  Table  I,  or  as  the  lower  value  if  a  range 
was  indicated, and  (ii)  only  if  with  that  choice  the  pin  will  serve  as  the 
softer  material  as  compared  to  the  average  or  upper  hardness  listed  for 
the  disk  material  in  Table  I.  When  even  so  the  disk  material  is  found  to 
be  softer  the  corresponding  box  is  left  blank.  Note  the  very  pronounced 
changes  of  flash  temperature  expected  on  interchanging  pin  and  disk,  by 
comparing  paired  values  along  diagional  lines  from  bottom  left  to  top  right.  For 
example,  g  ^  changes  from  13  for  a  magnesium  pin  sliding  on  a  copper  disk  to 
81  for  a  copper  pin  on  a  magnesium  disk.  Hence  the  occasionally  considerable 
changes  in  wear  rates  on  such  interchanges  become  understandable.  Signifi¬ 
cantly,  the  highest  value  of  all  arises  for  a  steel  pin  on  a  granite  (flint) 
substrate.  Thus  the  choice  of  materials  in  flintlock  rifles  was  most  apt 
and,  according  to  this  table,  perhaps  the  only  feasible  one 


SUBTRATE  SIDE 


PN  *  7rR2  p 

PT  *  flPU  =/17TR2p 


AT0/ATch  =  AT0/(^HS //olCl)  «  vr/(l/Z0-Xr) 


I 


Part  II  Fig.  3 


ATe/ATch=  ATe/(/i.Hs//D1c1)  -  vr/(l/Z0S  +  Xr/S 
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l’.i  rt  1 1  Fig .  4 


ATe/ATch  -  ATa/(/iHs//»lc1)  =  vr/(l/Z0S 


Part  II  Fig.  5 


ATe/ATch  =  ATe/t/iHg/^c,)  =  vr /(1/Z0S  +  Xr/S 
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Abstract 

A  simplified  model,  evolved  in  conjunction  with  previous  theoretical  work 
by  Blok  and  Jaeger,  yields  for  the  flash  temperature  at  elliptical  contact 
spots,  when  the  heat  is  uniformly  developed  over  the  contact  area, 
iTe  =  ( qR/x j)/( 1/ZqS  +  x r/ SQ ) ,  where  the  symbols  have  the  meaning  indicated  in 
the  list  of  symbols.  The  values  of  the  functions  ZQ  and  S  depend  on  relative 
velocity,  vf  =  v/vQ,  such  that  Zo(vr<2)  *  1/(1  +  vp/3)  ,  and  ZQ(vr>2) 

(9/8)/(/^  +  l/v8)  ,  and  S(vr<  1)  =  SQ  =  (4/3)//(  l+e3/4/3)  (l  +  l/3e3/4) ,  and 

S ( v  >  1 )  -  (e*  +  eV/8vr)/(  1  +  e^Z  //8vr) .  These  expressions  are  accurate  to 
within  a  few  percent,  and  they  are  represented  graphically  for  some  selected 
values  by  means  of  curves.  They  can  be  generalized  for  the  case  that  the 
contact  spot  moves  with  velocity  v^  and  v^  relative  to  the  two  sides, 

respectively,  as  iTg  =  q  R  /  [  '*  1/ZQ(vrl)S(vrl )  +  >  2/,Zo^vr2^S^vr2^  ’  where  vrl  = 

v ,/(<,/ R )  and  v  „  =  v0/(^-/R). 
i  i  r  c  c  c 


(i) 


Symbols 


a  radius  of  circular  contact  spot 

A  length  of  contact  spot  in  direction  of  motion 

B  breadth  of  contact  spot  at  right  angles  to  direction  of  motion 

c  specific  heat 

C  adjustable  factor  in  the  order  of  unity 

e  =  A/'B  axis  ratio  of  elliptical  contact  spot 

F  area  of  contact  spot 

K  adjustable  factor  in  the  order  of  unity 

M  adjustable  factor  in  the  order  of  unity 

q  heat  input  per  unit  area  and  time 

Q  rate  of  heat  input 

fWF/-r  characteristic  dimension  of  contact  spot 

S  shape  factor 

SQ  shape  factor  at  zero  velocity  (and  up  to  vr=l) 

t  time 

T  temperature 

iT  average  temperature  of  contact  spot  above  ambient 
v  velocity 

vQ  characteristic  velocity  (v  =  <^//F/tt) 

V  heated  volume 

Z  parameter  (function  of  contact  spot  shape  and  velocity) 

needed  for  computing  flash  temperatures . 

Z  for  circular  contact  spots 

<-\/oc  thermal  diffusivity 

heat  conductivity 

o  mechanical  density 


(ii) 


underlies  eq.  24,  is  significant  for  vr  =  1.  The  reason  for  this  breakdown  is 
that  forward  diffusion  (and  correspondingly ,  of  course,  back  diffusion)  was 
neglected.  Therefore  the  cooling  of  the  spots  is  increasingly  underestimated 
as  speeds  decrease.  Indeed,  also  comparison  with  Jaeger's  data  points  in  Fig. 
5  suggests  that  vr  =  1  is  already  mildly  below  the  applicability  of  the 
theory,  since  these  points  lie  noticeably  below  the  curve  for  e  <  1  and  above 
it  for  e  5  1. 

Therefore,  somewhat  arbitrarily,  for  the  present  purposes  the  best  values 
of  S  as  computed  from  eq.  24  were  judged  to  be  those  for  vr  =  9/8,  i.e. 

Sg/8  =  e1/4  ( l+l/3)/( l+e3/4/3)  =  e1/4(4/3)/( l+e3/4/3)  (29) 

Correspondingly ,  SQ  was  found  by  symmetrizing  Sg^g  as 

SQ  =  /S9/8(e)  x  Sg/8(l/e)  =  (4/3)/i+e3/4/3)  ( l+l/3e3/4)  (30a) 

The  function  S  has  been  entered  in  Fig.  5. 
o 

Expansion  to  the  Case  that  the  Contact  Spot  Moves  Relative  to  Both  Sides 

The  above  flash  temperatures  are  established  almost  immediately  at  the 
onset  of  motion.  Specifically,  after  the  spot  has  moved  once  through  its 
whole  length  the  equilibrium  flash  temperature  is  essentially  established 
(10),  as  would  indeed  be  expected  on  the  basis  of  the  theoretical  derivations 
presented  in  this  paper.  However,  it  is  possible,  and  even  probable  in  gears, 
that  the  spot  moves  relative  to  both  sides,  i.e.  at  velocity  Vj  in  regard  to 
material  1  and  v0  in  reaard  to  material  2. 

In  that  case  we  may  generalize  eq.  1  for  circular  contact  spots  to 

■To  =  qR/C'  1/Zo(vrl)  "  x2/Zo(vr2)] 
with  vrl  =  Vj/(«-  j/R)  and  vr2  -  v^/Uo/R). 
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(31a) 

(31b) 


Clearly,  as  given  in  eq.  28  is  an  excellent  representati on  of  ZQ  for 
v  <  2  as  will  be  seen  from  Fig.  4. 
b)  Elliptical  Spots 

The  remaining  task  is  to  treat  elliptical  spots  at  rest  or  very  low 
speeds,  and  to  this  end  we  must  first  obtain  a  shape  factor  for  the  stationary 
case  and  low  speeds.  This  is  necessary  not  only  for  cases  in  which  vr  <  2, 
but  because  in  regard  to  the  asperity  side  the  spot  is  stationary  at  all 
times,  and  we  must  know  in  which  manner  the  heat  flow  is  affected  by  shape 
changes . 

We  arrive  at  an  approximation  for  the  sought  value  S(v=o)  =  S0  by 
considering  Fig.  5  as  follows:  The  values  of  S  are  clearly  unsymmetrical 
about  e  =  1,  so  as  to  be  higher  for  e  than  1/e  if  e>  1 .  This  is  so  because 
spots  elongated  in  sliding  direction  leave  a  narrower  track  than  those 
elongated  at  right  angles  to  the  direction  of  motion,  and  thus  are  not  as 
effectively  cooled  by  the  substrate.  However,  both  lateral  and  forward  heat 
diffusion  increase  with  the  length  of  the  spot  perimeter  and  therefore,  at  low 
speeds,  elliptical  contact  spots  are  cooler  than  circular  ones  of  same  area, 
independent  of  axis  orientation.  Correspondingly ,  for  e  +  0  as  well  as  e  +  ■» 
we  find  S  tending  to  zero  at  all  speeds. 

The  approximation  of  S  according  to  eq.  24,  i.e.  on  the  basis  of  the 
high-speed  approximation,  evidently  is  still  quite  good  for  vr  =  1,  when  the 
maximum  value  of  S  (and  thus  maximum  temperature)  is  still  found  at  e  >  1,  as 
physically  it  must  at  all  speeds  v  >0.  However,  for  values  of  vf  below  1/2, eq. 
24  yields  maxima  of  S  for  e  values  smaller  than  1,  which  is  certainly  wrong. 
This  indication  of  a  noticeable  error  already  at  v^= 1  is  not  surprising  since 
we  know  from  Fig.  4  that  the  inaccuracy  of  the  high-speed  approximation,  which 


implicit  assumption  that  rounded  and  rectangular  spots  can  be  treated  as 
almost  the  same  is  vindicated  by  success,  since  S  does  not  deviate  much  from 
unity  for  most  reasonable  values  of  e. 

Refinement  of  the  Theory  for  v^  <2 

a)  Circular  Spots 

The  heat  diffusion  pattern  about  a  circular  spot  at  large  distances 
compared  to  a  is  essentially  spherical  at  all  speeds  on  the  asperity  side,  and 
is  so  on  the  substrate  side  provided  that  v  <<  v  .  Even  so,  the  interfacial 
temperature  will  decrease  with  speed  no  matter  how  small.  On  the  basis  of  the 
expected  distortion  of  the  diffusion  pattern  we  may  now  estimate  the  effect  of 
small  speeds  on  the  flash  temperature  as  follows:  On  the  substrate  side,  the 
effect  of  slow  movements  of  the  contact  spot  is  to  elongate  the  pattern  in  the 
direction  of  motion  linearly  with  the  velocity,  thereby  to  spread  out  the 
generated  heat  over  a  larger  volume,  and  thus  to  increase  the  heat  flow  into 
the  substrate  from  its  value  of  =  hA^aT  (see  eq.  2b).  The  thought  lies 
close  at  hand  to  represent  this  increase  in  by  the  factor  (1+Mvr),  where  M 
is  an  adjustable  factor  which  cannot  be  far  from  unity.  We  thus  obtain, 
instead  of  eqs.  2  to  4, 

Qx  =  irXjaaT  (1+Mvr)  (26) 

2 

and,  with  Q  =  =  na  q, 

O 

na  q  =  [  i  x^a(  1  +  Mvr)  +  irA^a]  aT  .  (27) 

Comparison  with  Fig.  1  yields  M  =  1/3  and  thus 

Z  Q  =  1/(1  +  Mvr)  =  1/(1  +  vr/3)  (28) 
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fitting  the  adjustable  constant  C  by  comparing  Z^e  with  1^Q  for  the  circle 
(i.e.  e  =  1)  at  high  speed.  Namely,  comparison  between  (21)  and  (22)  and  use 
of  eq.  8  renders 

Zhe  =  (9/8)ei/[(l+e3/*/C&Tr  )  /T  ]  =  (9/8)e*/(  ^+e3/4/C£  (23a) 

or  for  e  =  1  and  C  =  2 

Zh0  =  (9/8)/(^  +  1/vB  )  .  (23b) 

Fig.  4  indicates  the  result  for  V  for  the  circular  spot  according 
to  eq.  23b  for  our  initial  tentative  choice  of  C  =  2.  This  indeed  yields  a 
remarkably  good  fit.  Thus  encouraged  we  introduce  the  shape  factor 

S  =  Zhe/ZhQ  =  e^l+l/^/U+e374/^)  (24) 

so  that  we  may  write,  for  ellipsoidal  rather  than  circular  spots,  for 

vr>>(yx1)2  , 


SATo  . 


(25) 


Fig.  5  shows  the  shape  factor  for  a  range  of  vr  values,  including  vr  =  1, 
together  with  the  five  values  contained  in  Jaeger's  work  (10).  These  are 
remarkably  close  to  those  computed  from  eq.  24  for  vr  =  1.  This  is  presumably 
the  choice  made  by  Jaeger  since  he  did  not  specify  the  value  of  v^.  At  any 
rate  there  either  was  an  oversight  or  an  error  since  the  shape  f  ctor 
certainly  depends  on  speed.  Also,  Jaeger's  computations  of  S  (as  indeed  also 
for  Z)  were  made  for  rectangular  rather  than  round  spots,  and  in  this  much  our 
derived  values  cannot  be  entirely  accurate.  However,  as  seer,  from  Fig.  5,  our 
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above  estimate  of  the  average  temperature  in  the  quarter-ell ipsoid. 

Before  embarking  on  a  consideration  of  the  effect  of  sideways  diffusion, 
reflected  in  the  oB/B  term,  we  examine  the  case  of  very  high  speeds  in  which 
it  is  negligible.  From  eq.  (17)  this  is  found  as 

AThe  =  3Kqa/(4o1c16Z  v)  =  (3K/4)q/S/( /x^c^  /v) 

=  (3K//S)q  A  ek  /(/T^  ft  )  . 

AThe  =  (3K/,/5)q,/K  e^/Zx^jCj  kv  . 

Through  comparison  of  eq.  (18)  with  eq.  (5a)  we  deduce  that 
K  =  ( /8/3 ) ( 9/8 )  =  3//8  =  1.06  ,  and  that  in  the  high-speed  approximation  the 
flash  temperature  depends  on  ellipticity,  e,  as 

lThe'lTho  =  e  1/4  '  <>«> 

Returning,  then,  to  eq.  (17)  we  rewrite  it  as 

ATe/AThe  =  e1/4/(l+6B/B)  .  (20) 

Now  with  eq.  14b,  eq.  20  transforms  into 

AThe/AThQ  =  e*/(  l+e3/4/C  /^T)  .  (21) 

It  is  convenient  to  regard  the  above  dependence  of  the  flash  temperature  on 
the  parameter  C  as  due  to  a  shape  dependence  of  the  function  Z,  writing 

AThe  =  (22) 

in  parallelity  with  the  first  part  of  eq.  5a.  This  affords  an  opportunity  of 
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6  B  =  6  Z/C 


(14a) 


•  B/B  =  e3/4/(C/^T), 


(14b) 


expecting  C  to  be  in  the  order  of  2.  (We  shall  soon  see  that  C  =  2  is  an 
excellent  guess,  when  we  compare  the  present  approximation  with  Jaeger's 
results  (ref.  10)). 

In  the  present  approximation  we  neglect  forward  diffusion,  although  it 
does  occur  and  a  diffusion  front  precedes  the  spot  at  a  distance  of 

5 A  =  A/8/e  vr.  (15) 

The  reason  for  this  procedure  is  that  the  forward  diffused  heat  which  we 
neglect  at  the  front  end  of  the  ellipsoid,  is  balanced  by  the  heat  that  had 
already  been  deposited  ahead  of  the  midline,  into  the  quarter-ellipsoid,  by 
the  spot  in  its  initial  position.  However,  this  approximation ,  as  well  as 
that  for  <5 B  according  to  eq.  14,  and  indeed  the  entire  model,  becomes  invalid 
at  speeds  below  vo,  or  so. 

The  volume  of  the  heated  quarter-el  1 ipsoid  under  consideration  (Fig.  3) 


V  =  (tt/6 )AB ( 1  +<5 B/B )<5 z  =  ( 2/3 ) F ( 1+6B/B)  6z 


whence,  by  the  use  of  eq.  12 


>Te  =  *KQin/plClVe  =  3KqA/[2p1c1(l+cB/B)6z  v] 


where  K  is  an  adjustable  factor  near  unity.  The  factor  K  is  necessary  since 
one  can  hardly  expect  the  interfacial  temperature  to  be  exactly  equal  to  the 
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presented  by  him,  was  abandoned  after  many  hours  and  after  the  University 
Computer  Center,  having  had  the  project  under  advisement  for  several  days, 
declined  the  task.  Therefore,  we  shall  give  a  simple  estimate,  applicable  to 
vr>  2  or  perhaps  even  vr  >  1,  by  refining  the  model  already  used  to  derive 
eq.  5: 

By  geometry,  about  one  half  of  the  heat 

Qin  =  AFq/v,  (12) 

developed  while  the  contact  slides  once  through  its  whole  length  A,  is 
distributed  in  a  volume  whose  shape  is  included  in  a  diffusion  front  as 
sketched  in  Fig.  3.  The  average  temperature  in  the  outlined  quarter-ellipsoid 
is  representative  of  the  contact  spot  temperature,  within  some  factor  not  far 
from  unity,  according  to  the  previous  derivation  of  ATho 

(eq.  5).  The  sideways  extension  of  the  ellipsoid  beyond  the  width  of  the 
track  swept  out,  5B,  is  due  to  lateral  diffusion,  and  its  short  axis,  6z,  is 
due  to  inward  diffusion.  Much  as  before,  seeing  that  inward  diffusion  occurs 
during  time  interval  A/v,  we  estimate  6z  as 

<5z  =  /k^A/v  =  A/(e*/2vr  )  .  (13) 

Clearly,  6B  is  some  fraction  of  6z  since,  firstly,  lateral  diffusion  of 
significance  begins  only  after  the  contact  spot  has  reached  the  mid-way 
position  and,  secondly,  the  heat  content  outside  the  geometrical  limits  of  the 
swept  track  of  width  B  is  fed  only  very  ineffectively  after  the  perimeter  of 
the  spot  shrinks  away  from  them  once  again.  On  the  other  hand,  6B  consists  of 
two  strips,  one  at  either  side  of  the  track. 


Further  we  redefine  the  characteristic  velocity  v  as 

J  o 

vQ  =  vF7  it  =  <j/R  (10) 

in  general,  which  we  met  already  for  circular  spots,  in  particular,  as 
vo  ■  *l/a- 

It  is  mostly  helpful  and  instructive  to  use  relative  velocities 

v  =  v/v  ,  rather  than  absolute  values.  The  significance  of  v  is  that,  in  the 
r  o  3  o 

substrate,  at  velocity  vQ,  the  contact  spot  and  diffusion  front  emitted  at 
time  t=0  both  traverse  the  characteristic  distance,  R,  in  the  characteristic 


time  interval 


St  =  R/vq  =  R  /  <p 


v  St  =  St  =  R. 
o  1 


(11a) 


(lib) 


Essentially,  therefore,  at  speeds  much  below  vQ  heat  diffusion  in  the  sub¬ 
strate  keeps  up  with  the  contact  spot.  Conversely,  at  v  >>  vq,  as  far  as  the 
contact  spot  is  concerned,  the  heat  is  left  behind  in  the  substrate  where  it 
was  first  deposited.  It  is  thus  vq  in  regard  to  which  v  has  to  be  either 
small  or  large  to  generate  the  low  and  high  speed  cases  already  discussed. 


The  Effect  of  Contact  Spot  Shape  at  v  >2  When  the  Asperity  Side  is  an 


Insulator 

No  straight-forward ,  simple  treatment  of  the  effect  of  contact  spot  shape 
on  flash  temperature  is  available.  Jaeger  (10),  Holm  (4)  and  Shobert  (5)  have 
considered  it  mathematically  to  some  extent,  but  their  equations  are  cumber¬ 
some,  to  say  the  least,  and  it  is  difficult  to  find  specific  answers  to  actual 
cases  from  these.  Indeed,  an  initial  attempt  at  computing  a  number  of  re¬ 
presentative  values  from  the  equations  given  by  Jaeger  (10),  beyond  the  few 
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is  deposited.  Underneath  the  spot,  the  heat  meanwhile  had  time  to  diffuse 


only  through  a  thickness  in  the  order  of  /<^2a/v  <<  a,  and  thus  it  is  essen¬ 
tially  dispersed  in  a  volume 

V  =  4a2  /kj  2a/v  .  (7) 

The  interfacial  temperature  rise  must  be  comparable  to  the  average 
temperature  rise  calculated  for  that  volume  since,  on  the  one  hand,  not  all  of 
the  heat  energy  is  confined  to  the  volume  V  and,  on  the  other  hand,  the 
interface  is  hotter  than  the  average  temperature  in  V.  Therefore 

ATho  *  Qin  /(piciV)  =  qvi/(/x1P1c1  /v  )  (5b) 

which  is  to  be  compared  with  eq.  5a.  We  thus  found  in  the  high-speed  limit 
for  circular  contacts  Zho  .1  /T~7av  which  is  fortuitously  close  to  the 
actual  answer  embodied  in  Fig.  2,  namely 

ZhQ  =  (9/8)  =  (9/8 )/^JVQ  =  (9/8)/^  .  (8) 

Characteristic  Dimension,  Velocity,  and  Distance 
Between  Spot  and  Thermal  Diffusion  Front 

In  general  a  contact  spot  will  not  be  circular  but  roughly  elliptical, 
with  axes  lengths  A  and  B,  say,  in  the  direction  parallel  and  at  right  angles 
to  the  direction  of  motion.  The  ellipticity  of  the  spot  ^ ha  11  be  called 
e  =  A/B  and  its  area  F,  so  that  we  may  define  the  characteristic  dimension  of 


we  find 


aTq  =  qa/(Aj+  \2)  (4) 

for  the  contact  spot  at  rest.  Comparing  eqs.  1  and  4  shows  thus  that  at  rest, 
or  at  small  speeds,  ZQ  =  1,  in  agreement  with  Fig  2. 

c)  The  Limiting  Case  of  Very  High  Speeds 

When  the  asperity  slides  over  the  substrate,  it  continuously  encounters 
cool  substrate  and  leaves  a  trail  of  warmed  material  behind,  much  like  the 
contrail  of  an  airplane.  At  very  low  speeds  the  extra  heat  energy  needed  to 
heat  up  that  trail  is  negligibly  small  compared  to  the  radial  heat  diffusion 
of  the  stationary  case  just  discussed.  But  as  the  heated  volume  increases,  so 
does  the  heat  left  behind  in  the  trail.  Thus  with  rising  velocity  the  contact 
spot  is  increasingly  cooled  by  the  substrate,  until  the  heat  absorbed  in  the 
trail  becomes  dominant  at  high  speeds.  At  given  q,  therefore,  AT  decreases 
with  speed,  mediated  through  the  corresponding  decrease  of  Z  and,  depending  on 
A^/Ag.  sooner  or  later  the  heat  flow  into  the  asperity  side  becomes  negligi¬ 
ble. 

For  that  case,  i.e.  Q2  < (  Q^,  from  eq.  1  in  conjunction  with  Fig.  2  and 
introducing  the  relative  speed  vr  =  v/vQ  =  v/(<^/a),  one  finds  directly 

AT  (vr>>l)  =  AThQ  =  qaZhQ/X1  =  (9/8)  qa/A^vy  =  (9/8)  q/a/(  A^Cj  *V).  (5a) 

We  may  readily  derive  this  result  semi-quantitatively  as  follows: 

Consider  the  heat  flow  into  the  substrate  while  the  contact  spot  moves  at  high 

o 

speed  through  the  distance  2a,  in  time  2 a/v,  generating  (2a)  of  new  trail 
area,  while  the  heat 


4 


where  the  subscripts  1  and  2  designate  a  substrate  (1)  and  flat  asperity  (2), 
respectively,  as  sketched  in  Fig.  1.  Here  x  stands  for  heat  conductivity,  q 
is  the  amount  of  heat  developed  per  unit  area  of  contact  spot  per  unit  time, 
and  ZQ  is  a  numerically  known  function  of  velocity.  The  subscript  o  is 
carried  as  a  reminder  that  the  function  ZQ  and  the  above  temperature  rise  Tq 
apply  to  circular  contact  spots  only. 

Eq.  1  was  derived  from  the  early  work  of  Block  (8,9)  as  elaborated  by 
Jaeger  (10,11).  The  function  ZQ  is  indicated  in  Fig  2  as  extracted  from  Fig. 
7  of  ref.  10. 

b)  The  Limiting  Case  of  Very  Low  Speeds 

The  asperity  side  is  stationary  in  regard  to  the  contact  spot  at  all 
times.  If  it  is  circular  of  radius  a,  heat  diffuses  into  the  material  of  the 
asperity  side  at  the  rate 

Q2  =  itX^T,  (2a) 

if  the  interface  is  at  temperature  a T  above  ambient,  regardless  of  relative 
velocity.  To  the  extent  that  the  bulk  temperatures  of  the  two  sides  are 
equal,  as  is  assumed  throughout,  and  that  ideal  thermal  contact  exists  between 
them  at  the  contact  spot,  then  at  rest,  i.e.  for  velocity  v  =  0,  the  heat  flow 
into  the  substrate  must  be  similarly 

=  nX  ^ad  T.  ( 2b ) 

Since  the  total  heat  flow  is  the  sum  of  those  two  parts,  i.e. 

Qj  +  Q2  =  Q  =  *a\T  (x  ^x  2)  =  qra2,  (3) 
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to  a  lack  of  simple  methods  to  compute,  or  at  least  estimate,  them.  In  regard 
to  the  latter,  one  is  painfully  aware  that  the  scientist  will  ordinarily 
disdain  numerical  estimates  and  regard  quantitative  data  the  more  highly,  the 
more  accurately  they  are  known.  In  regard  to  flash  temperatures  that  would  be 
an  unjustifiable  attitude,  though.  In  actual  fact,  numbers  and  average  sizes 
of  contact  spots  change  from  moment  to  moment,  and  any  individual  contact 
spot,  as  long  or  briefly  as  it  may  last,  does  not  only  change  its  size  and 
shape,  and  hence  its  temperature,  within  short  time  intervals,  but  also 
experiences  erratic  changes  in  the  momentary  temperature  distribution  over  its 
area.  Therefore,  while  it  is  possible  to  calculate  the  temperature  distribu¬ 
tion  over  any  specific  contact  spot  once  its  size,  speed,  shape  and  the  nature 
and  thickness  of  intervening  surface  layers  or  lubrication  layers  are  known, 
what  would  be  the  benefit  thereof  unless  one  felt  confident  that  one  dealt 
with  the  representative  contact  spot?  It  is  the  average  contact  spot, 
therefore,  after  which  one  should  enquire,  and  one  must  realize  that  actual 
flash  temperatures  will  widely  fluctuate  about  the  value  calculated  for  it. 

In  this  spirit,  then,  we  shall  develop  some  simple  equations  to  aid  in 
such  calculations.  In  Parts  II  and  III  we  shall  then  consider  practical 
applications  and  give  some  sample  values,  expanding  on  previous  work  (6-7) 

Basic  Equations 

a)  The  General  Equation  for  Circular  Contact  Spots 

For  circular  contact  spots  of  radius  a,  one  may  write  for  the  average 


flash  temperature  (6,7) 
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Introduction 

The  temperature  rise  at  contact  spots  during  sliding  between  two  solids 
with  or  without  simultaneous  current  flow  has  great  significance  in  a  number 
of  ways,  and  some  well-known  papers  addressed  to  the  non-mathematician  are 
available  (e.g.  refs.  1-3).  Thus  temperature  spikes  at  contact  spots  induce 
thermal  stresses  which  on  repeated  contacts  can  cause  fatiguing;  they  mostly 
cause  superficial  softening  and  thus  reduction  of  the  friction;  they  can  cause 
the  formation  or  dissolution  of  precipitates;  can  anneal  out  surface  modifica¬ 
tions  such  as  through  ion  implantation;  will  generally  enhance  surface  re¬ 
actions,  especially  oxidation;  will  enhance  diffusion  into  and  out  of  surfaces 
so  as  to  cause  surface  decarburization  of  steels,  for  example;  and  under 
suitable  circumstances  can  cause  phase  transformations.  An  important  factor 
in  the  enumerated  temperature  effects,  especially  in  regard  to  chemical 
reactions,  diffusion,  and  changes  in  precipitates,  is  the  fact  that  cooling  is 
typically  so  rapid  that  the  thermal  equilibrium  vacancy  concentration  of  the 
flash  temperature  is  essentially  frozen  in.  Correspondingly  these  phenomena 
can  be  enhanced  much  beyond  the  expectations  on  the  basis  of  the  time-temper¬ 
ature  history  at  the  contact  spots  alone. 

In  spite  of  the  above,  and  in  spite  of  the  already  cited  papers  plus  the 
important  work  by  Holm  (4)  and  Shobert  (5),  much  uncertainty  and  con  fusion 
exists  regarding  the  magnitude  of  flash  temperatures.  Presumably  this  is  due 
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SB  increase  of  track  width  at  trailing  end  of  contact  spot  due  to 
lateral  diffusion  of  heat 

Sz  thickness  of  heated  layer  underneath  contact  spot 
Subscripts 

e  pertaining  to  elliptic  contact  spots 

h  pertaining  to  high  speeds  (v  >>  1) 

l  pertaining  to  low  speeds  and  to  rest  (v  <(  1) 

o  pertaining  to  circular  contact  spots  (except  in  the  case  of  v 

which  applies  also  to  elliptical  spots)  0 

1  pertaining  to  the  substrate 

2  pertaining  to  the  asperity  side 


Further,  on  the  basis  of  the  considerations  in  the  preceding  section  we  obtain 
for  elliptical  contact  spots  in  that  case 


ATe  =  qR/[x1/Z0(vrl)S(vrl)  +  yZ0(vr2)S(vr2)].  (32) 

Summarizing  Conclusions 

The  accuracy  with  which  the  present  simple  theory  reproduces  Jaeger' s( 10) 
theoretical  results  is  most  reassuring  and  inspires  confidence  that  the  model 
effectively  represents  the  physics  involved.  If  we  accept  the  flash 
temperature  of  the  circular  spot  at  rest  with  insulating  asperity-side,  of 
value  q //F/tt/x,  ,  as  basic,  it  is  reasonable  and  convenient  to  express  all 
flash  temperatures  as  fractions  thereof.  Namely,  under  no  circumstances  does 
the  flash  temperature  rise  above  the  mark  of  qa/x^  =  q/ZF/ir/X^. 

For  circular  spots  in  motion  with  a  non-insulating  asperity  side  of  thermal 
conductivity  x2,  naming  X^/X^  =  Xr,  we  may  then  write 

AT0/(q/F7 V  /Xj)  =  l/(l/Z0+xr).  (33) 

As  was  shown,  in  accordance  with  eqs.  23  b  and  28  and  demonstrated  in  Fig.  4, 
Z  is  well  represented  by 

Zo(vr<2)  =  Z^  =  1/(1  +  vr/3)  (28) 


and 


Zo(V2)  =  Zho  =  9/8/(^V  +  1Z/§  )  •  (23c) 

The  maximum  error  with  this  approximation  is  about  3%  at  vr  =  2,  as  seen  in 
Fig.  4. 

Fig.  6  shows  aT  for  a  few  values  of  x  between  0  (i.e.  insulating 
a  o  r  3 

asperity)  and  xr  =  10  (i.e.  the  thermal  conductivity  of  the  substrate  is  only 
10%  of  that  of  the  asperity  side).  From  this  in  conjunction  with  Fig.  2  it  is 


evident  that  the  high-speed  limit,  in  which  the  bulk  of  the  heat  is  absorbed 
by  the  substrate,  is  already  fairly  well  approximated  at  vr  =•  2  for  =  0, 
but  only  for  much  higher  speeds  at  larger  values  of  a.  In  fact  for  the 

O 

high-speed  approximation  to  be  useable,  it  must  be  v  >JA  . 

The  effect  of  ellipticity  is  expressed  in  t<  nns  of  the  shape  factor  S. 

For  vr  <  1  it  is  approximately  equal  to  SQ,  as  given  in  eq.  30a  or,  in  a 
slightly  different  but  not  superior  formulation, 

S  ( v f <1 )  =  SQ  =  (4/3) A  10/9)  +  e3/4/3  +  l/3e3/4  .  (30b) 

Thus,  at  slow  speeds  the  shape  factor  applying  to  the  substrate  is  S  ,  the 
same  as  for  the  asperity  side.  For  elliptical  shapes  (including,  of  course, 
circles  with  e  =  1)  at  low  speed  we  therefore  obtain 

ATfe/(q*r7^  /Aj)  -  S0/(l/ZJlo+Ar)  =  SQ/(1  +  v/3  +  ^  .  (34) 

For  higher  speeds  though,  as  certainly  for  vr  >  1,  the  shape  factor  S  applies 
to  the  substrate  but  SQ  continues  to  apply  to  the  asperity  side,  so  that 

ATe/(qvF7^  /Aj)  -  1/(1/Z0S  +  Ar/S0).  (35) 

Figs.  7,  8  and  9  show  the  flash  temperature  computed  according  to  eqs.  34 
and  35,  whereby  between  v„  =  0.2  and  v  =  1  the  transition  between  S  and  S 
has  been  smoothed.  It  is  interesting  to  note  how  e  values  above  unity  lift 
the  flash  temperatures  above  that  for  circular  spots  at  higher  speeds,  and  how 
those  speeds  at  which  ATg  is  larger  than  aTq  increase  with  increasing  e 
values  and,  for  same  e,  with  increasing  A  Conversely,  at  low  speeds  any 
ellipticity  reduces  the  flash  temperature,  and  does  so  by  the  same  amount  for 
e  and  1/e.  For e  <  1  the  factor  by  which  the  temperature  is  reduced  increases 
with  increasing  speed.  Even  so,  the  effect  of  ellipticity,  whether  e  is  above  or 


below  unity,  is  non-major  for  all  speeds  below  vr~  100,  staying  within  a 
factor  of  two,  more  or  less,  of  that  for  circular  spots,  at  all  values  of 
vr  below  100  for  0.1  <  e  <  10. 
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Legends 

Fig.l)  Mathematical  model  of  a  contact  spot,  consisting  of  a  fiat  circular  or 
elliptical  asperity  of  material  2  with  thermal  conductivity  ^ 
sliding  with  velocity  v  on  a  semi-infinite  substrate  of  material  1 
with  thermal  conductivity 

Fig. 2)  The  parameter  ZQ  as  a  function  of  relative  velocity  vr  =  v/vQ  for  a 

circular  contact  spot  as  extracted  from  Jaeger's  calculation  according 
to  Fig.  7  of  ref.  10.  Also  indicated  are  the  high-speed  limit 
Zh0  =  (9/8)//v^  and  the  low  speed  limit  of  unity  given  in  ref.  10. 

Fig. 3)  Model  of  the  volume  (roughly  a  quarter-ellipsoid)  about  an  elliptical 
contact  spot  (with  axes  A  and  B)  which  receives  about  one  half  of  the 
heat  energy  developed  while  the  spot  slides  once  through  its  length  A. 

Fig. 4)  The  function  ZhQ  =  (9/8)/(/vy  +  1//8  ),  applicable  for  vr  s  2,  and  the 
function  Z^  =  l/(l+vr/3)  applicable  for  vr<  2,  in  accordance  with  eqs. 
23b  and  28  . 

Fig.  5)  The  shape  factor  S  =  e4  ( l+l//8v  )/( l+e^//8v^) ,  and  the  shape  factor 

SQ  =  4/3//l+e3/4/3)( l+l/3e3//4)  for  stationary  spots  according  to 
eqs.  24  and  30a,  respectively.  SQ  is  applicable  to  vr  <  1  and  S  to 
vr  >  1.  Included  are  values  given  by  Jaeger  (ref.  10)  for  rectangular 
contact  spots  (symbol  ■  )presumed  to  be  applicable  to  vr  =  1. 

The  shape  factor  corrects  the  function  ZQ  for  application  to  elliptical 
contact  spots.  Namely  for  vr>l  it  is  Z^g  -  SZhQ  and  for  vr  £  1  it  is 

Z„  =  S  Z  .where  Z.  and  In  are  as  shown  in  Fig.  4. 
ste  oo  ho  ^o 
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Fig. 6)  The  flash  temperature  of  a  circular  contact  spot  in  units  of  q/F /Tr/ A^ , 
i.e.  as  fraction  of  the  temperature  which  that  spot  would  have  at  rest 
when  the  asperity  side  is  an  insulator,  as  a  function  of  relative 
velocity  for  different  values  of  A^  =  A^/A^.  Large  values  of  A^  mean 
that  the  asperity  side  is  a  much  better  conductor  than  the  substrate 
and  thus  conducts  away  much  of  the  heat. 

Fig. 7)  The  flash  temperature  for  contact  spots  of  different  shapes  in  units 
of  q/F/r  /\  =  qR/A^  (i.e.  as  fraction  of  the  temperature  which  a 
circular  spot  of  same  area  would  have  at  rest  when  the  asperity  side 
is  an  insulator)  for  A^  =  A^/A.  =  0  (meaning  in  this  case  the  asperity 
side  j_s  an  insulator).  El  1 iptici ties  e  =  A/B  >  1  mean  spots  elongated 
in  sliding  direction  and  e  <  1  spots  flattened  in  sliding  direction. 

Fig. 8)  As  Fig.  7  but  for  A ^  =  1,  meaning  that  the  heat  conductivities  on  the 
two  sides  are  equal . 

Fig.S)  As  Fig.  7  but  for  Ar  =  10,  i.e.  the  case  in  which  the  asperity  side  is 
a  ten  times  better  conductor  than  the  substrate. 
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ABSTRACT 

Approximate  solutions,  reliable  within  a  few  per¬ 
cent,  have  been  obtained  for  the  theoretical  flash 
temperatures  at  circular  and  elliptical  contact  spots 
between  two  homogeneous  naterials  when  heat  is  evolved 
at  the  (assuned  mathsnatical  ly  flat  and  unoontaminated) 
interface  at  a  uniform  rate,  either  independent  of 
velocity  as  in  the  case  of  Joule  heat,  or  in  proportion 
with  the  relative  speed,  he.  as  is  the  case  for  frict¬ 
ion  heat.  The  parameters  appearing  in  the  corresponding 
equations  have  been  expressed  in  terms  of  known  mater¬ 
ials  properties  plus  the  number  of  contact  spots  and 
their  axis  ratio.  The  accuracy  with  which  the  flash 
temperatures  can  be  known  in  practice  is  limited  by 
uncertainty  regarding  those  parameters,  and  the  ideal¬ 
izations  embodied  in  the  model.  Those  theoretical  un¬ 
certainties  are  probably  no  larger  than  the  actual 
momentary  fluctuations  of  the  conditions  at  the  contact 
spots. 


a  Characteristic  dimension  of  current-carrying 

circular  contact  spots,  i.e.  a-spots 

gA  Tbtal  load-bearing  area  in  the  elastic  case 

p/V=P/Hs  Tbtal  Load  Bearing  area  in  the  plastic  case 

A  Length  of  load-bearing  rontact  spot  in  direct¬ 

ion  of  motion 

B  Breadth  of  load-bearing  contact  spot  at  right 

angles  to  A 
c  Specific  heat 

d  Mechanical  density 

e=A/B  Axis  ratio  of  elliptical  load-bearing  contact 
spot 


E  Young's  modulus 

f  (Z,S)  Describesthevelocity  andshape  dependence  of 
the  flash  temperature:  ATe=  AT0f(Z,S) 

F=h(AB/4  Area  of  load-bearing  contact  spot 
Hs  Mechanical  hardness  of  the  softer  side 

I  Tbtal  current  flowing  between  the  two  sides 

N  Number  of  contact  spots 

P  "Load",  i.e.  normal  force  between  the  two  sides 

q  Heat  input  per  unit  time  and  contact  spot  area 

Q  Rate  of  total  heat  generation 

Characteristic  dimension  of  load-bearing  con¬ 
tact  spot 

t2  Radius  of  curvature  of  asperity  making  elastic 

contact  with  the  substrate 
R  Electrical  resistance 

S  Shape  factor 

S0  Shape  factor  at  zero  velocity  and  up  to  vr=l 
T  Temperature 

T  Average  temperature  of  contact  spot  above  ambient 

T0=q r/A^  Reference  temperature  (circular  spot  at  v=0 
when  the  asperity  is  thermally  insulating) 

v  Sliding  velocity  of  the  asperity  (side  2)  over 

the  substrate  (side  1) 

v0  Characteristic  velocity  (Vo'Xj/r  or  v^Kj/a, 

respectively) 

vr=v/v0  Relative  velocity  between  the  two  sides 

Z  Function  needed  for  computing  the  velocity 

dependence  of  the  flash  temperature 

Z0  Z  for  circular  contact  spots 


* 


T 


T 


K=A/dc  Thermal  diffusivltiy 
Heat  oondutivity 

Heat  conductivity  of  the  asperity  side  re¬ 
lative  to  the  substrate 
Coefficient  of  friction 
■>  Poisson's  ratio 

^  Electrical  resistivity  of  the  sliding  materials 

Cf  Film  resistivity  (  electrical  resistance  of  unit 
area  of  film) 

Subscripts 

a  Pertaining  to  a-spots 

c  Pertaining  to  constriction  resistance 

e  Pertaining  to  elliptical  contact  spots 

E  Pertaining  to  elastically  deformedcontact  spots 

f  Pertaining  to  electrical  film  resistance 

F  Pertaining  to  friction 

J  Pertaining  to  Joule  heat 

P  Pertaining  to  plastically  deformed  contact  spots 

r  Standing  for  "relative" 

T  Pertaining  to  current  conductionviaelectron 

tunneling  through  surface  films 

1  Pertaining  to  side  1  of  tiie  sliding  pair,  the 

substrate 

2  Pertaining  to  side  2  of  the  sliding  pair,  the 

asperity 


FUNDAMENTAL  RELATIONSHIPS 

In  a  previous  set  of  two  papers  HI  the  following 
expression  was  derived  for  the  flash  temperature  at  an 
elliptical  contact  spot  between  two  semi- inf  ini  te 
materials  as  shown  in  Fig.  1,  based  on  the  fundamental 
earlier  work  fcy  Blok  ill  and  jaegor  <3.41. 

5  qf/iV[zolv[)s(vt)l  +Vso}  (1) 

where 

vr  =  v/v0  =  v/{jq/r)  (2) 

with  subscripts  1  and  2  referring  to  the  two  sides  as 
indicated  in  Fig.l  and  in  the  list  of  symbols,  where 
also  the  other  symbols  are  explained. 

The  function  ZQ  was  numerically  kncwn  from  ref.3, 
but  no  analytical  expression  for  it,  except  that  for 
very  low  speeds  it  becomes  unity,  and  at  very  high 
speeds  approximates  (9/8) /•/v^. 

By  means  of  a  simple  model,  in  which  the  diffusion 
distances  of  the  heat  energy  are  considered, the  fol¬ 
lowing  approximations  were  derived  and  found  to  compare 
to  Jaeger's  numerical  values  within  a  few  percent  over 
the  whole  range  of  speeds  ill: 


Z0(vt42)  =  V(1  +  Vj/3) 

(3a) 

and 

Z0(vr22)  -  (9/8)/(>|Tc  +  VVF )  . 

(3b) 

Using  the  sane  approach  the  shape  factor,  S,  was 
found  to  be  approximately  given  by 

S(vril)  =  S0  =  (4/3)/[ (1  +  e3/4/3)(l  +  l/3e3/4]^2  (4a) 
and 

S(vr21)  =  (e1/4  +  eV4/V£rrVU  +  e3/4/V8vq).  (4b) 

In  order  to  emplcy  egs.l  to  4  above  for  the  evalu¬ 
ation  of  flash  temperatures  in  actual  cases,  one  must 
knew  the  various  parameters  appearing  therein.  Except 
for  q  and  r,  these  are  obvious  and  will  ordinarily  be 
known  within  reasonably  narrow  limits.  The  discussion 
therefore  now  turns  to  the  determination  of  q  and  r.  In 
this,  one  may  treat  different  additive  contributions  to 
the  total  heat  evolution  separately,  at  least  to  a 
first-order  approximation,  since  the  excess  of  the 
flash  temperature  above  ambient  is  simply  proportional 
to  q,  as  seen.  Correspondingly,  we  may  write 

<3  =  <3f  +  Sj  (5) 

where  qF  and  qj  signify  the  heat  evolution  due  to 
friction  and  joule  heat,  respectively.  However,  since 
q  is  the  rate  of  heat  evolution  £e_i  unit  area  _oj£  con¬ 
tact  spot,  one  must  generally  not  only  know  how  much 
total  heat  is  evolved  in  order  to  know  q,  but  also  the 
size  of  the  spots,  i.e.  r  or  a.  As  we  6hall  see,  this 
causes  some  complications;. 


EVALUATION  CF  C  AM)  a.  THE  CHARACTERISTIC  SPOT  SIZES 

(i)  Hie  Dimension  rr  of  Plastic  Contact  Spots 

When  the  load-bearing  contact  spots  ate  plastical¬ 
ly  defotmed,  the  average  normal  pressure  over  the  whole 
area  of  real  contact,  independent  of  the  number  and 
shape  of  participating  contact  spots,  approximates  Hs, 
the  (Meyer)  harchess  of  the  softer  of  the  two  sides,  be 
it  the  asperity  or  substrate  side. 

In  that  case,  then,  the  total  load  between  the  two 
sides,  P,  is  supported  by  the  total  area  pA  of  the 
plastically  deformed  load-bearing  contact  spots,  say,  N 
in  nunber,  so  that 

pA  =  P/Hs  =  Nn  pr2  (6a) 

or 

pr  =  VP/iNHs  .  (6b) 


coincide  with  load-bearing  areas.  Overwhelmingly  (but 
as  we  shall  see  not  exclusively)  current  passes  only 
where  load-bearing  contact  is  established.  Yet  on 
account  of  uneven  aoverage  by  surface  films  of  varied 
nature,  enormous  local  fluctuations  of  contact  resist¬ 
ance  and,  hence,  of  current  density  can  exist,  such  as 
at  cracks  in  oxide  leyers.  Correspondingly,  no  blanket 
statements  as  to  the  values  of  j  can  be  made,  but  each 
case  must  be  considered  individually.  Only  in  the 
presence  of  very  uniform  surface  coverage  with  either  a 
conducting  film  or  a  "tunneling  film",  Le.  a  film  thin 
enough  to  permit  electron  tunneling,  can  general  equat¬ 
ions  be  applied. 

If  either  type  of  uniform  film  is  present,  and  the 
oontact  spots  are  plastically  deformed,  the  associated 
value  of  the  a-spot  radius  is  sinply 

pa  =  pr  ,  (8) 


(ii)  ihe  Dimension  Ec_  of  Elastic  Contact  Spots 

In  the  majority  of  cases,  when  wear  occurs,  the 
contact  spots  are  plastic,  in  accordance  with  the  above 
derivation.  However,  at  very  light  loads,  or  when  two 
very  hard  or  brittle  materials  contact,  or  else  if 
repeated  contact  is  made  at  the  same  spots,  the  behav¬ 
ior  can  be  elastic.  In  that  case,  following  R.  Holm  (5) 
and  again  for  N  contact  spots  under  total  load  P, 

Er  =  { (3P/4N)  ( (1-  vf )/Ej+  (1-M|)/E2J  r2)1/3  (7a) 

where  -0  stands  for  Poisson's  ratio,  E  for  Young's  modu¬ 
lus,  and  r2  is  the  radius  of  curvature  of  the  asperity. 

Cn  the  assumption  that  neither  nor  *72  differ 
too  much  from  1/3,  and  with  l/f^  +  l/Pj  «  1/Pgve 
this  somewhat  cumbersome  expression  becomes 

Er  =  (2Pr2/3NEave)1/3  ,  (7b) 

and  correspondingly  the  total  load-bearing  area  in  the 
elastic  case  is  approximately  given  by 

gA  =  Nirgr2  =  2.4  (NP^/E av|)1/3  .  (7c) 

Strictly  speaking,  the  model  of  Fig.l  on  which 
eqs.  1  to  4  are  based  is  incompatible  with  an  el  as  tic 
contact  spot  caused  by  a  spherical  asperity  meeting  a 
half-plane,  described  by  the  Hertzian  eq.7.  However, 
the  heat  flow  between  the  two  sides  is  negligibly 
affected  if  Fig.2  is  substituted  for  Fig.  1,  as  implied 
by  eq.7.  The  error  resulting  therefrom  is  certainly 
smaller  than  that  due  to  any  of  the  other  causes, 
unless  r2  is  comparable  in  magnitude  to  A. 

Illil-Ihe  .Blmensiong_of  Current  Carrying  Contact  Spots 

The  areas  of  current-carrying  contact  spots,  gen¬ 
erally  referred  to  as  "a-spots"  J5L  do  not  necessarily 


and  similarly  for  comparatively  large  elastic  contact 
spots 

E3  =  Er  ■ 

Finally,  it  is  possible  that  the  contact  spots  are 
so  small  relative  to  the  radius  of  curvature  of  the 
asperity  (i.e.  for  Ao<r2  in  Fig.2)  that  the  tunneling 
current  about  the  periphery  of  the  load-bearing  spots 
becomes  significant.  Namely,  electrons  tunnel  not  only 
through  very  thin  films  but  similarly  through  any  nar¬ 
row  gap  between  two  materials,  and  thus  also  through 
the  annular  wedge-diaped  gap  about  the  periphery  of  the 
contact  spots.  This  effect  may  be  enhanced  by  elastic 
tv- r h -n  through  adhesive  forces  acting  through  the 
Same  gap,  leading  to  a  slight  expansion  of  the  con¬ 
tacting  area.  As  a  tesult,  current  will  flew  through  a 
somewhat  larger  area  than  necessary  for  elastic  load- 
bearing,  but  the  corresponding  increase  of  a  above  is 
negligible,  indeed  undiscoverable,  except  for  extranely 
snail  elastic  contact  spots. 

In  previous  papers  evidence  for  peripheral  tunnel¬ 
ing  has  been  discovered  and  the  corresponding  theoreti¬ 
cal  considerations  have  been  presented  (6-10).  The 
appropriate  value  of  j  is  151 

Ta  =  (Er2  +  2sr2)1/2  (10) 

with  £  the  gap  width  (before  deformation  through  atjie- 
sive  forces,  if  any)  through  which  tunneling  takes 
place.  In  the  cases  studied  in  refs.  £.  to  fi,  the 
paraneter  £  was  found  to  be  near  5  angstroms. 

E¥AUiAXIQtLQELs rr  DUE  TO  FRICTIONAL  HEAT 

When  two  materials  elide  past  each  other  with 
velocity  v  under  the  load  P  with  a  friction  coefficient 
y  ,  heat  is  developed  at  the  rate 


Qp  =  ^PV  .  (11) 

Since  this  heat  evolution  takes  place  at  the  contact 
spots  we  find,  for  the  subscript  i  standing  for  either 
Por  E,  i.e.  the  plastic  or  elastic  case, 


Similarly,  when  with  entirely  uniform  surface 
coverage  the  contact  spots  are  elastic,  so  that  eq.  7 
may  be  used,  one  obtains  approximately 

Efb  Er  =  0-153  I2  <Sf E^e/PNtj  .  (20) 


pqp  =  Qf/N-fpr2 


iqF  pr  =  Pv/Nrr  jr  . 

For  the  plastic  case,  thus 
pqF  Pr  =  ^vifPHg/Nti  , 


and  for  the  elastic  case 


e3f  Er  =  0.364 ^v(P2Eave/r2N2) 1/3 


At  or  very  near  the  interface  Joule  heat  is  devel¬ 
oped  from  two  causes.  One  of  these  is  the  electrical 
resistance  of  the  ubiquitous  surface  films,  and  be  they 
only  monontolecular  absorption  layera  If  these  are  not 
present,  the  two  sides  weld  together,  and  the  theory 
breaks  down  in  any  event  on  account  of  "galling"  or 
"seizing". 

Designating  the  specific  film  resistance  by  the 
symbol  6f,  the  total  Joule  heat  generated  ie 


The  second  contribution  to  Joule  heat  at  or  neat 
the  interface  is  due  to  the  constriction  resistance. 
This  arises  because  all  of  the  current  has  to  pass 
through  the  a- spots  and  thus  the  current  flow  lines  are 
strongly  oonstricted  at  the  interface,  the  effect  is 
much  as  if  the  two  solids  were  connected  by  Na  wires  of 
the  same  diameter  as  the  a-spots  and  a  length  compar¬ 
able  to  their  radius,  but  unlike  the  film  resistance  it 
depends  on  the  shape  of  the  a-spots.  Approximately,  the 
constriction  resistance  of  Na  elliptical  a-spots  with 
axis  ratio  e  -  A/B  and  area  IT  ^a2  each  is,  according 
to  eqs.4.23  and  4.27  of  ref.5.# 

i*c  =  (fi+?2>Un(2^e)/Ve]/Na-|Tia  (21a) 

where  once  again  the  subscript  i  can  signify  P,  E  or  T, 

while  ^ is  the  electrical  resistivity  of  the  sliding 
materials.  Hcwever,  for  eil.4 

i«c  -  <?1+  ?2>/4Na",ia  (2U,) 

isebetter  approximation  which  indeed  is  exact  for  e=l. 
Correspondingly,  it  is 

icQj  =  *2i*c  (22> 


fQj  =  l2Pf  «  l2df/Na-<  1a2  (16) 

where  ja  can  stand  for  pa,  ga  or  .ja,  as  the  case  may 
be,  and  the  subscript  a  at  Na  shall  draw  attention  to 
the  fact  that  the  number  of  a-spots  can  be  larger  or 
smaller  than  the  nunber  of  load-bearing  spots,  and  this 
by  factors  as  large  as  ten  either  way,  and  perhaps 
larger  yet. 

Correspondingly , , 

fqj=  I^f/IHa^a2)2  d?) 


f^J  ia 


For  the  particular  case  that  a  uniform  film  exists 
overthewhole  surface,  we  have  Na=N  with  pa  =pr  for 
plastic  contact  spots  and  ga  ■  gt  for  elastic  spots  in 
accordance  with  equations  6a  and  h,  If  so,  we  have 

Pf'i;  Pr  ”  I2  <5f(Hg  /P'Nn)1/2  .  (19) 


irfj  s  icQ/NaTfi*2 


icflj  ia  =  I2(fr+?2><l*(2^>//e]/<Na-t'ia>2  .  <24> 

For  the  case  of  uniform  surface  coverage  with  a 
conducting  film  or  a  tunneling  film  we  thus  obtain,  for 
plastic  contact  spots, 

Pdlj  Pr  =  l2<?l  +?2>tln(24e)/Ve)(Hs/PUN)  .  (25) 

Similarly,  for  elastic  contact  6pots  obtain 

g^qj  Er  ■=  0.133  I2('5,1+V2)Iln(2ve)/;e)(Eave/N2Pr2)2/3 

(26) 

We  shall  not  specifically  discuss  the  case  of 
reduction  of  joule  heat  through  peripheral  tunneling 
since  this  is  of  very  specialized  interest  only,  is 
ordinarily  a  minor  effect,  and  has  been  discussed  in 
refs.  8. to  1.0.  Therefore,  in  summary  of  this  section, 
then,  the  heat  flow  density  due  to  electric  currents  is 

i'lj  "  ifTj  +  icflj  ■ 


(19) 


We  m^  rewrite  eq.l  as 

Xre=Kr0t(Z,S)  (27a) 

where 

AT0=  qr/Aj  «  q/F/n  '/Aj  (27b) 

and 

f(Z,S)  =  (l/Z0(vt)S(vt)  +  Vso)_1*  <27c) 

In  this  general  formulation,  ATe  is  the  sum  of  all 
contributions  to  the  flash  temperature  which  simply 
superimpose  (at  least  to  a  first  approximation)  as  q  is 
the  sun  of  all  contributions.  Thus 

We  =  ATeF+  ATgj  (27d) 

and 

q  =  qF  +  qj  =  qF  +  (27e) 

as  discussed  in  the  previous  section,  while 

At  =  V^l  (27f) 

and  vt=v/vQ=v/(  Xj/r)  in  accordance  with  eq.2.  As 
before,  r  is  the  characteristic  dimension  of  the  con¬ 
tact  spot,  nanely 


were  given  above  in  eqs.6-10. 

As  we  saw,  it  is  not  generally  possible  to  give 
values  for  qj  except  when  the  surface  film  has  uniform 
electrical  resistance  everywhere.  For  that  particular 
case,  based  on  eqs.13,  18  and  24  obtain 

ATe  =  [f(Z,S)/Aj]  [^Pv/(Nirir)  +  I2 Gf/tN2  K2  ir3) 

+I2(J>1  +  f2Hln(2Ve)//i]/(NTfir)2j  (28) 

where  is  either  pr  or  Er. 

Hence  f ran  eqs  14,  19  and  25 

p£Te  =  [f(Z,S)/,l1]^WFHs/Nti'+  I^fH^/P^N)1/2 

+  i2<9i+  92>  Un(2V5)/V5)  (HB/PT1N)  j  (29) 

and  from  eqs  15,  20  and  26 

gATe  «  [f(Z,S)/A1){0.364p(P2Eave/r2N2)1/3 

+  0.153  l2  (3fEave/Pm2 

+  0.133  I2(f1+92)(ln(2Ve)/7e](Eave/N2Pr2)2/3i[  (30) 
while 

pvr  *=  vVP/NTTHe  /Ki  (31) 


r  =  Vf/  U  , 


(27g) 


and 


of  which  the  specific  values  pr,  Er,  p a,  and  .ja 


Evr  =  v(2Pr2/3NEove)1/3/K1  . 


(32) 


FIG. 3 


DISCUSSION 

T0*  qr/Aj  is  the  flash  temperature  of  a  circ¬ 
ular  contact  spot  at  rest,  comprising  a  thermally  insu¬ 
lating  asperity  side,  but  having  the  same  area,  located 
on  the  same  substrate,  and  experiencing  the  same  rate 
of  heat  input,  as  the  contact  spot  mder  oonsideratioa 
The  function  f(Z,S)  corrects  this  basic  flash  temper¬ 
ature  in  accordance  with  the  actually  existing  modific¬ 
ations  due  to  the  facts  that  the  actual  oontact  spot  is 
generally  (i)  in  motion,  (ii)  elliptical,  and  (iii) 
that  the  asperity  side  is  ordinarily  not  a  thermal 
insulator  so  that  not  all  of  the  heat  flows  into  the 
substrate  (although  the  fraction  that  flows  in  that 
direction  increases  steadily  with  increasing  speed 
because  more  and  more  of  the  heat  is  expended  in  warm¬ 
ing  up  a  surface  layer  of  the  fresh  track  on  the  sub¬ 
strate)  . 

All  three  of  the  named  effects  diminish  the  flash 
temperature,  so  that  A.T0  is  the  highest  flash 
temperature  that  can  be  generated  at  a  spot  of  the  area 
F  ■  n  r^  situated  on  the  given  substrate  materiaijand 
heated  ty  heat  input  q  per  uiit  area-  Correspondingly, 
the  function  f(Z,S),  shewn  for  At=  0  and  for  \  *  10 
in  Figs. 3  and  4,  is  neverjlarger  than  unity,  and 
6teadily  decreases  with  rising  speed  and  heat 
conductivity  of  the  asperity  side. 

In  regard  to  ellipticity  f(Z,S)  is  not  symmetrical 
with  respect  to  e  =  A/B.  This  is  for  the  reason  that 
a  oontact  spot  stretched  out  in  the  direction  of  motion 
leaves  tehind  a  rarrewer  warmed  track  on  the  substrate 


than  a  spot  of  same  area  which  is  wide  at  right  angles 
to  the  direction  of  motion;  and  thus  the  flash  temper¬ 
ature  of  the  former  is  higher  than  of  the  latter. 

While  evidently  none  of  the  derivations  or 
expressions  given  above  are  inherently  difficult,  their 
practical  application  is  still  cumbersome,  albeit  ty 
far  more  managable  than  in  the  few  previously  available 
treatments.  Moreover,  the  functional  dependence  of  the 
flash  temperature  on  velocity  was  knewn  only  imperfect¬ 
ly  and  its  shape  dependence  even  less,  although  one 
should  hasten  to  add  that  Holm  (5)  has  by  a  different 
route  arrived  at  another,  and  perhaps  more  accurate, 
method  of  evaluating  flash  temperatures  under  a  wide 
range  of  circimstances. 

The  variations  in  which  the  equations  can  be  writ¬ 
ten  are  manifold,  and  in  semewhat  different  formulat¬ 
ions  they  have  been  assembled  in  form  of  a  table  in 
refs.  £  and  id,  where  also  seme  nunerical  examples  for 
the  flash  temperature  are  given  for  the  case  of  circul¬ 
ar  a-spots. 

Cue  has  no  reason  to  be  surprised  that  the  equat¬ 
ions  are  so  involved:  The  nunber  of  physical  parameters 
and  variables  which  in  actuality  determine  the  flash 
tenperature  is  considerable,  and  somehow  the  equations 
have  to  accomodate  that  fact. 

We  have  already  pointed  out  at  the  beginning,  that 
the  reliability  of  the  obtained  numerical  values  can  be 
no  better  than  the  input  parameters  plus  the 
approximations  inherent  in  the  model.  Yet,  results 


that  have  variously  been  computed  tor  particular  cases 
(8.101  give  the  impression  of  being  gratlfyingly  close 
to  fact,  wherein  as  a  rule  of  thumb  N  =  10  is  a  good 
choice  unlesr  one  knows  the  number  of  contact  spots 
more  reliably  from  some  other  source. 
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